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Abstract 

We show that an ordinary semiconducting thin film with spin-orbit coupling can, under ap- 
propriate circumstances, be in a quantum topologically ordered state supporting exotic Majorana 
excitations which follow non-Abelian statistics. The key to the quantum topological order is 
the coexistence of spin-orbit coupling with proximity-induced s-wave superconductivity and an 
externally-induced Zeeman coupling of the spins. For the Zeeman coupling below a critical value, 
the system is a non-topological (proximity-induced) s-wave superconductor. However, for a range 
of Zeeman coupling above the critical value, the lowest energy excited state inside a vortex is a 
zero-energy Majorana fermion state. The system, thus, has entered into a non-Abclian s-wave 
superconducting state via a topological quantum phase transition (TQPT) tuned by the Zeeman 
coupling. In the topological phase, since the time reversal symmetry is explicitly broken by the 
Zeeman term in the Hamiltonian, the edge of the film constitutes a chiral Majorana wire. Just like 
the s-wave superconductivity, the Zeeman coupling can also be proximity-induced in the film by an 
adjacent magnetic insulator. Wc show this by an explicit model tight-binding calculation for both 
types of proximity effects in the heterostructure geometry. Here we show that the same TQPT 
can be accessed by varying the interface transparency between the film and the superconductor. 
For the transparency below (above) a critical value, the system is a topological (regular) s-wave 
superconductor. In the one-dimensional version of the same structure and for the Zeeman cou- 
pling above the critical value, there are localized Majorana zero-energy modes at the two ends of a 
semiconducting quantum nanowire. In this case, the Zeeman coupling can be induced more easily 
by an external magnetic field parallel to the wire, obviating the need for a magnetic insulator. We 
show that, despite the fact that the superconducting pair potential in the nanowire is explicitly 
s-wave, tunneling of electrons to the ends of the wire reveals a pronounced zero-bias peak. Such 
a peak is absent when the Zeeman coupling is below its critical value, i.e., the nanowire is in 
the non-topological s-wave superconducting state. We argue that the observation of this zero-bias 
tunneling peak in the semiconductor nanowire is possibly the simplest and clearest experiment 
proposed so far to unambiguously detect a Majorana fermion mode in a condensed matter system. 

PACS numbers: 03.67.Lx, 71.10.Pm, 74.45.-Fc 



2 



I. INTRODUCTION 



Particle statistics of a collection of indistinguishable particles is a genuinely quantum 
mechanical concept without any classical analog. In spatial dimensions three and above, 
pairwise interchange of particle coordinates in a many-body system is equivalent to a sim- 
ple permutation of the coordinates. Consequently, each interchange has the effect of either 
a change of sign (fermion) or no change at all (boson) on the many-body quantum wave 
function. In (2 -|- 1) dimensions, however, exchanges and permutations are not necessarily 
equivalent. In this case, under simple interchange of the particle coordinates, the cor- 
responding space-time trajectories can form non-trivial braids in the (2 + 1) -dimensional 
space-time. ^ Consequently, in (2 -|- 1) dimensions, particles can have quantum statistics 
strikingly different from the statistics of bosons and fermions. 

A straightforward extension of the statistics of bosons and fermions is the Abelian anyonic 
statistics, in which the many body wave-function, under pairwise exchange of the particle co- 
ordinates, picks up a phase 6, which can take any value between (bosons) and vr (fermions). 
Since a phase factor is only a one-dimensional representation of the braid group in 2D, the 
statistics is still Abelian. On the other hand, if the many body ground state wave function 
happens to be a linear combination of states from a degenerate subspace, a pairwise ex- 
change of the particles can unitarily rotate the wave function in the ground state subspace. 
In this case, the effect of exchanging the particle positions is an operation on the wave func- 
tion vector by a unitary matrix representing this rotation. Consequently, the statistics is 
non- Abelian,^ and the corresponding system is a non- Abelian quantum system. It has been 
proposed that such systems, if the ground state subspace is concurrently protected by an 
energy gap, can be used as a fault-tolerant platform for topological quantum computation 
(TQC).HHS! 

One important class of non- Abelian quantum systems, sometimes referred to as the Ising 
topological class,'' is characterized by quasiparticle excitations called Majorana fermions, 
which involve no energy cost (when the mutual separation among the excitations is large). 
The second quantized operators, 7^, corresponding to these zero energy excitations are self- 
hermitian, 7J = 7j, which is in striking contrast to ordinary fermionic (or bosonic) operators 
for which Cj 7^ cj. However, since 7^ and 7j anticommute when i ^ j, they retain some prop- 
erties of ordinary fermion operators as well. The Majorana fermions, which are actually more 
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like half-fermions, were envisionecP by E. Majorana in 1935 as fundamental constituents of 
nature [e.g. neutrinos are thought to be Majorana, rather than Dirac, fermions). Majorana 
modes are intriguing because each Majorana particle is its own anti-particle unlike Dirac 
fermions where electrons and positrons (or holes) are distinct. Although the emergence 
of Majorana excitations, which are effectively fractionalized objects (anyons) obeying non- 
Abelian anyonic statistics, in solid state systems is by itself an extraordinary phenomenon, 
a great deal of attention has also been focused on them because of the possibility of carrying 
out fault tolerant TQC in two dimensional systems using these Majorana particles. TQC, 
in contrast to ordinary quantum computation, would not require any quantum error correc- 
tion since the Majorana excitations are immune to local noise by virtue of their non-local 
topological nature.l^ The direct experimental observation of Majorana modes in solid state 
systems would therefore be a remarkable breakthrough both from the perspective of funda- 
mental physics of fractional statistics in nature and the technological perspective of building 
a working quantum computer. It is therefore not surprising that there has been recent 
resurgence of immense interest for the experimental realization (and detection) of Majorana 
fermions in solid state systems. Recently, some exotic ordered states in condensed matter 
systems, such as the Pfaffian states in fractional quantum Hall (FQH) systems,'^^'^ p-wave 
superconductors/superfiuids,'^^'!^^ theoretical models that can be potentially simulated in 
cold atom optical lattice systems,'^^^ as well as the surface state of a topological insulator 
(TI) or related systemd^^HlSI have been discussed as systems which can support Majorana 
fermions as the lowest energy excitations. In the context of optical lattice systems, it has 
also been proposed that a 2D + ipy superfluid can be realized using only s-wave Fesh- 
bach resonance modified by the topological Berry phases arising from artificially generated 
spin-orbit coupling.lSSl 

It has been shown recently ^ that even a regular semiconducting film with a sizable 
Rashba-type spin-orbit coupling, such as InGaAs thin films, can host, under suitable con- 
ditions, Majorana fermions as low energy excitations. Since the basic effects behind the 
emergence of the Majorana fermion excitations - spin-orbit coupling, s-wave superconduc- 
tivity, and Zeeman splitting - are physically well-understood and experimentally known to 
occur in many solid state materials, the proposed semiconductor heterostructure is possi- 
bly one of the simplest condensed matter systems supporting Majorana quasiparticles and 
non-Abelian quantum order. By an analysis of the real-space Bogoliubov-de Gennes (BdG) 
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equations for a vortex in the semiconductor thin film, in which s-wave superconductivity 
and a Zeeman sphtting are proximity-induced (Fig. (la)), it has been shown that the lowest 
energy quasiparticle excitation in the vortex core is a zero-energy Majorana fermion mode. 
This real space analysis has also been supported by a momentum space analysis in the form 
of an index theoremP^l analogous to such a treatment in the context of one-dimensional Dirac 
theory.'SsEll Here a comment about the various means to induce a Zeeman splitting in the 
semiconductor thin film is in order. Note that when the spin-orbit coupling is of the Rashba 
type, we require a Zeeman splitting which is perpendicular to the plane of the film (Zeeman 
splitting parallel to the film does not produce a gap in the one-electron band-structure, a firm 
requirement of our non-Abelian state) .'^^'^ Inducing such a splitting by applying a strong 
perpendicular magnetic field is not convenient, because the magnetic field will give rise to 
unwanted order parameter defects such as vortices. It is for this reason that we propose 
to induce the Zeeman splitting by the exchange proximity effect of an adjacent magnetic 
insulator (we ignore the small coupling of the spins in the semiconductor with the actual 
magnetic field of the magnetic insulator). More recently, it has been shown that, when the 
spin-orbit coupling also has a component which is of the Dresselhaus type, the appropriate 
Zeeman splitting can also be induced by applying an in-plane magnetic field.'^The Majorana 
mode is separated by a finite energy gap (so-called mini-gap) from the other conventional 
fermionic excited states in the vortex core. Thus, for a collection of well-separated vortices, 
the resulting degenerate ground state subspace is protected from the environment by the 
mini-gap. This enables the potential use of the semiconductor heterostructure in Fig. (la) 
in TQC. 

One of the main goals of the present paper is to provide the important mathematical 
details relevant to our solutions of the BdG equations in the semiconductor heterostructure. 
These mathematical details, which were left out in Ref. [51], are given in Sec. [II] and 
Sec. [Ill] below. It is important to note that, unlike the case of the surface of a 3D strong 
topological insulator adjacent to an s-wave superconductor,!^ the BdG equations in the 
spin-orbit coupled semiconductor are not exactly solvable. We therefore only show that, 
in a specified region of the parameter space, a single non-degenerate solution of the BdG 
equations, which is spatially localized around the vortex core, is allowed. We also show 
that the second quantized operator corresponding to such a solution is indeed a Majorana 
fermion operator. (In a subsequent section (Sec. [VI]) we confirm the existence of such zero- 
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energy Majorana fermion states localized at the vortex cores by a full numerical solution 
of the BdG equations set up on a sphere.) In the next few sections (Sees. [IV-IX]) we 
provide a comprehensive discussion of the interesting physics of non-Abelian topological 
order arising via the complex interplay of spin-orbit coupling, Zeeman splitting, and s-wave 
superconductivity externally induced in a host system. We also deduce the parameter space 
needed for the establishment of the non-Abelian order, as well as the associated TQPT 
accessed by tuning the system in or out of this parameter space. In Sec. [X] we study the 
superconducting and magnetic proximity effects in the host semiconductor thin film by a 
microscopic model tight-binding calculation. In the last part of the paper (Sec. [XI]), we 
consider a one-dimensional version of our proposed structure- a semiconducting nanowire 
with proximity induced s-wave superconductivity and a Zeeman splitting. We emphasize 
that the Zeeman splitting can now be induced by a magnetic field parallel to the length of 
the wire, because such a field does produce a gap in the one-electron band structure without 
producing unwanted excitations in the adjacent superconductor. This obviates the need for 
a nearby magnetic insulator. For the Zeeman splitting above a critical value, the wire is 
in a non-Abelian topological phase with zero-energy Majorana excitations at the ends. We 
propose a scanning tunneling experiment from the ends of the semiconducting nanowire as 
possibly the most realistic experiment proposed so far to detect a Majorana fermion state 
in a condensed matter system. 

We use the terminology Majorana particle or Majorana fermion or Majorana state or 
Majorana excitation or Majorana mode interchangeably in this article, all of them meaning 
precisely the same entity, namely, the non-degenerate zero-energy eigenstate {i.e. a solution 
of the BdG equations) at the vortex core of a spin-less chiral p-wave or other such topological 
superconductor. We emphasize that this object obeys the non-Abelian braiding statistics 
rather than ordinary fermionic statistics, and the Majorana particle is its own anti-particle 
in contrast to the ordinary Dirac fermions where electrons and holes (positrons) are distinct 
particle- hole conjugates of each other. A part of the results presented in this paper - 
approximate solutions of the BdG equations in the semiconductor heterostructure - has 
been published elsewhere. ^In Sees. [II, III] we provide all the mathematical details relevant 
to the solutions of the BdG equation which were left out in Ref. |5T]. Most of the results 
contained in the subsequent sections are new. Some additional mathematical details related 
to the solution of the BdG equations are relegated to the appendix. 
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FIG. 1: (a): The proposed heterostructure of a semiconductor (SM) sandwiched between an s-wave 
superconductor (SC) and a magnetic insulator (MI). In this geometry, the semiconducting film can 
support non-Abelian topological order, (b): Single-particle band-structure in the semiconducting 
film with and without the Zeeman splitting induced by the MI. Without the Zeeman splitting, the 
two spin-orbit shifted bands touch at = ky = \k\ = (red lines). Then, for any value of the 
chemical potential, the system has two Fermi surfaces. With a finite Zeeman splitting, the bands 
have an energy gap near kx = ky = \k\ = (blue lines). If the chemical potential lies in the gap, 
the system just has one Fermi surface (indicated by the dotted circle). 

II. HAMILTONIAN. 

The single-particle effective Hamiltonian Hq for the conduction band of a spin-orbit cou- 
pled semiconductor in contact with a magnetic insulator is given by (we set = 1 henceforth) 

Ho = ——-iJ.+V^a^+a{axp}-z. (1) 
2m* 

Here, m*, Vz and /i are the conduction-band effective mass of an electron, effective Zeeman 
coupling induced by proximity to a magnetic insulator (we neglect the direct coupling of 
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the electrons with the magnetic field from the magnetic insulator), and chemical potential, 
respectively. The coefficient a describes the strength of the Rashba spin-orbit coupling and 
(Tq, are the Pauli matrices. 

The proximity-induced superconductivity in the semiconductor can be described by the 
Hamiltonian, 

4 = j dv {A(r)4(r)c|(r) + H.c}, (2) 

where cj.(r) are the creation operators for electrons with spin a and A(r) is the proximity- 
induced gap. The pairing term Hp and the non-interacting part Hq can be combined to 
obtain the BCS mean-field Hamiltonian Hbcs = Hq + Hp. The excitation spectrum of this 
Hamiltonian is defined in terms of the Bogoliubov quasiparticle operators 



/*E",(r)4(r)..v(r)Ur) (3) 



which satisfy 

[Hbcs,1^]=EjI (4) 

Such a quasiparticle operator 7 can be used to construct excited states 7^|\l/o) with energy 
E + Eq from the ground state |^o) with energy Eq. The ground state |\E'o) is defined as 
the lowest energy state of the BCS Hamiltonian satisfying 7|\I^o) = 0. The equation for the 
quasiparticle operator, Eq. |4} can be re-written as the BdG equations in the Nambu basis, 

^0 A(r) \ 

\ ^(r) = E^{r). (5) 
A*(r) -cryH*ayJ 

Here, \l'(r) is the wave function in the Nambu spinor basis, \I'(r) = 

(n-f-(r), M4^(r), t>4^(r), — ff-(r))-^. Introducing the Pauli matrices Tq, in the Nambu space 

the Hamiltonian on the left hand side in Eq. (|5| can be written as 

2 

HBdG = [^-fi+V,a, + aiaxp)-z]T, + [A(r)r+ + h.c] (6) 
where r+ = rl = Ii±^. 



III. BDG EQUATIONS FOR A VORTEX. 

The single-particle Hamiltonian Hq can be written in polar coordinates as 
Ho = r]p'^ - fx + V^a^ + a{a x p) ■ z 

= --qV"^ - fi + V^az + t^icr+P- -(T-P+) (7) 
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where 77 = (7+ = ctL = o^^iOy and p+ = Px + iPy = e'^{-idr + ^de) and p_ = - zpj^ = 
e~^^{—idr — ^de). The full BdG Hamiltonian for an n-fold vortex can be written conveniently 
in the Nambu space as 

HbcLG = (-?7V^ - h)t^ + T^cr^ + ^-((7+p_ - (7_p+)t^ + A(r)[cos {ne)Tx + sin (n^)T2^]. (8) 

In order to diagonalize the above Hamiltonian it is convenient to note that the BdG 
Hamiltonian has a combined spin-orbit-pseudospin rotational symmetry. This symmetry 
can be expressed compactly by noting that HBdG commutes with the operator 



(9) 



Therefore, the eigenspinors of the BdG Hamiltonian can be taken to be eigenstates with 
eigenvalue = nij of the form 



(10) 



The above equation can be used to eliminate the angular degree of freedom 9 from the 
BdG equations as follows: 



HBdG,mj^mj{r) = -Emj^mjW- 



(11) 
(12) 



Here HBdG,mj = e-<"'^-'^^/2+"^-/2)e^^^^g*(mj-a./2+nr./2)e g independent. More specifically 



HBdG,mj ^ -{r]{d^ + -dr + 



1 (2mj - cr^ + TiTzf 



4^2 



) + fJ'}rz + V^a^ 



- - <y-}rzdr - 2 ^' 2 Hn^x- (13) 

Under the action of the particle-hole transformation operator, S = ayTyK, the mj spinor 
eigenstate with energy E transforms into a —mj eigenstate with energy —E because 

Therefore, a necessary condition for a non-degenerate, E — 0, Majorana state solution is 
that nij — 0. Prom here onwards we will write HBdG,mj=o — HBdG ^-nd ^^^=o(^) = ^(r). 



Single- valuedness of the spinor wave-functions in Eq. ( 10 ) requires that (n — 1)/2 must be an 
integer. Therefore, only odd vortices can have non-degenerate Majorana eigenstates. From 
here onwards, for the simplicity of discussion, we will consider only zero energy solutions at 
the cores of single-flux-quantum vortices (n = 1) . 

The BdG matrix HbcLG may be reduced to a real Hamiltonian by applying the cr^ rota- 
tion U = e^"'''^!^ as lisdc U^^HbcigU. The solutions of the resulting E = BdG equation 
HsdG^i''^) — must come in complex conjugate pairs \[^(r) and \l/*(r). Therefore the so- 
lutions \E'(r) can be required to be real without loss of generality. For such real solutions, 
it follows from the particle-hole symmetry of the BdG equations that ayTy'^ijr) is also a 
solution. Thus, any non-degenerate E = solution must be real and satisfy the property 
ayTy^ir) = A\E'(r). Moreover, because {cTyTyY = 1, the possible values of A are A = ±1. 

Using the relation = iXcTyTz, which follows from ayTy = A, the BdG Hamiltonian for a 
given value of A is of the form 



(y. 01. i~ -\- \ T — 1 

- 2 + ^-}^^^- - ^^^+^^ + "^-^2"^^^ ^ i\(yyT,/\{r). (14) 

The Hamiltonian in this limit does not couple the = ±1 sectors (electron and hole). 
This allows one to write the BdG differential equation in terms of only the electron-sector 
{tz = +1) of the spinor ^o('") = ■ '^^^ corresponding reduced BdG equations 

for a single vortex {n = 1) take the form of a 2 x 2 matrix differential equation: 



(^{d'r + -rdr) + Vz-f^ XA{r) + a{dr + h) \ 



^o(r)=0. (15) 



\^ -XA{r)-adr -^{dl + ldr-^)-Vz-ii] 

The hole part of the spinor is not independent and is constrained by the value of A such 
that v^{r) = Xu^{r) and v^^{r) = Xu^{r) and the Majorana spinor has the form \l/(r) = 
i^oir),tay^oirrf- 

We now approximate the radial dependence of A(r) by A(r) = for r < R and A(r) = A 
for r > R where R is of the order of the radius of a vortex core. In view of the topological 
stability of the putative Majorana zero-energy solution to local changes in the Hamiltonian, 
^ such an approximation can be made without any loss of generality. 
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A. Solution inside the vortex core. 



Inside the vortex core (r < R), which is the non-superconducting region (A(r) = 0), it 
is possible to construct exphcit analytic solutions \E'(r, z) to these equations in terms of the 
Bessel functions Jo{z) and Ji{z) as, 

*(r..)=r^'°'"'' |. (16) 

UiJi{zr) 



By substituting Eq. (16) into Eq. (15) we find that {u^,ui) and z satisfy 



{-r]z^ + V; - fi)u^Jo{zr) + zau^Jo^zr) 
zau^Ji{zr) + {—rjz^ ~ — ix)u^Ji{zr) 



0, (17) 



which implies 

-7]z'^ + Vz — ^J. za 

az —rjz'^ — Vz~fi 

Existence of solutions in terms of {u^, u^f requires that z satisfies the characteristic equation 




0. (18) 



of the matrix in Eq. (18) , which is given by 

{7]z^ -nf -V^ - z^c? = 0. (19) 

This is a quadratic equation in z'^ with real roots. Therefore, the solutions for z are either 
purely real or purely imaginary and come in pairs with opposite signs. A real root z = k 
of this equation corresponds to a crossing of some band at the Fermi level. For purely real 
roots z, only the solution \E'(r, z) corresponding to positive z is normalizable at the origin 
and therefore physically acceptable. On the other hand, the purely imaginary roots z = ±tk 
lead to a single real solution \E'(r, ik) = \E'(r, —ik). Thus we can see that in general there are 
2 linearly independent solutions. If all 4 solutions of z are real then these correspond to the 
2 Fermi surfaces obtained from the intersection of the bands with the Fermi level. If only 1 
pair of solutions is real then the imaginary pair corresponds to a decaying state. 

The BdG equations describing the proximity-induced supercondutivity at a TI/SC inter- 
face^ follow from the BdG equations for the present system, Eq. (15), by taking rj = 0. In 



11 



this case the matrix equation reduces to 

zn \ I lu \ 

= (20) 




This equation only has 1 pair of real solutions z = ii^/i^ — V^/a and therefore has only 1 
linearly independent solution in the core of the vortex. 

B. Solution outside the vortex core. 

The solution outside the vortex does not have a simple analytic form as the solution 
inside. Motivated by the large r asymptotic expansion for Bessel functions, for r > R we 
can consider a series expansion of the form 





(21) 

where p{x) are analytic power series in x. We expect to be able to close such a series of 



equations since the matrix in Eq. (15) only has derivatives and powers of 1/r. 

= (22) 



r 



1/2 




r]{-df. - Idr) + V,-^i AA + a{dr + D 

-AA - a{dr) vi-d'r - }dr + j,)-V,-fi 

^(-^r - 4^ + 2^9, - + v; - /i AA + a{dr + i - ^) W Pt{l/r) 

-AA - a{dr -i-z) i]{-dl + J, + 2zdr -z^)-V,-i2 [ p;(l/r) 



0. 



(23) 



As shown in the appendix, the last equation has a simple solution as a power-series in 
1/r which can be determined numerically. Moreover in this power series expansion we can 
determine the equation for the 0-th order term by formally setting 1/r = as below 



— r/^^ + \4 — AA — za 
— AA + za —riz"^ — Vz — ^ 




0. (24) 



Setting z = ik one can see that the matrix appearing in the above equation is related to 
the one determining the quasiparticle bandstructure from the BdG equations. The values 
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FIG. 2: Complex roots z = z/j + zz/ of the characteristic equation Eq. [25| shown in the complex 
plane. In the non-Abelian phase (Co = (A^ + — V^) < 0), 3 roots with positive real parts and 1 
root with negative real part for A = — 1 while there are only 2 roots on either side of the imaginary 
axis for Co > for A = — 1. The roots in the A = 1 channel are the negative of the roots in the 
A = 1 channel. 

of z consistent with the above equation are determined by setting 



Det 



-rfz^ + Vz — jJL AA — za 
— AA + za —riz^ — K — A* 

= {r]z' - - + {za T A)^ = 0. (25) 
The two famihes of solution for A = ±1 are related simply by flipping the sign of z. The 



sign of the product of the roots Zn of Eq. 25 is given by S* = sgn (Ilnl-^n)) = sgn (Co) where 
Cq = [fi"^ + A^ — yj) is the polynomial evaluated at 2; = 0. The parity of the number of 
normalizable solutions with Re{zn) > is given by P = sgn (fl^ _Re(z„)). Since Eq. 25 is 
real, complex roots z„ occur in conjugate pairs. Therefore, complex roots cannot affect the 
sign of either 5* or P. It follows that S = P. 



Therefore, the condition \ Vz\ > y A^ + /i^ implies that there are an odd number of roots 
with positive real parts. Specifically, as seen in Fig. (|2|), if Co < 0, there are 3 roots on one 
side of the imaginary axis and 1 root on the other side. Similarly for Cq > we must have 
2 solutions on each side of the imaginary axis. A slightly different version of this argument 
has previously been presented.!^ 
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C. Matching Boundary conditions at the edge of the vortex core. 



As discussed before, for + /i^ > V|, one of the channels labeled by ± has a solution 
with 3 decaying (negative real parts) solutions and 1 growing solution. Out of these 3 
decaying solutions, one is purely real and negative and the other 2 are complex conjugate 



with negative real parts. An observation that can be made by considering Eq. (23) is 
that if (p-f-(l/r), p4^(l/r))^ corresponds to a value z then (p;^(l/r), p^(l/r))'^ corresponds to 
an eigenvalue z* . Thus from one pair of complex conjugate decaying eignevalues we can 
construct a pair of real solutions 



Wt,{l,2)('^) 



S(i,2){e" 





(26) 



where Si = 1 and 8^ = i with the solution 1 corresponding to + and 2 to — . 
The non-degenerate real eigenvalue already corresponds to a real eigenvector 

pt(l/r) 

On the other hand, for r < R we expect a 2 parameter family with the general solution 
given by 




(27) 




+ as 



Matching the gradient and the wave-function at r = i? we get 





(28) 



/ MR) \ 

ViiR) 

drV^{R) 



04 



f Wt,4(^) \ 
\ drUl^4{R) j 



I ^tAR) \ 



+ ^5 



f n^AR) \ 
UijiR) 

drU^j{R) 
\drU^j{R) J 



(29) 



nAR) 

drU^^5{R) 
\drU^,5{R) J 

Together with the normalization constraint on the global wave-function, this leads to 5 
equations in 5 variables, which leads to a unique solution for the Majorana mode in the 
case Co < 0. However, for the other case with Co > 0, there are only 2 decaying modes 
outside the vortex core. The existence of a Majorana mode would then require us to satisfy 
5 equations with 4 variables. Such a problem in general is over constrained and no Majorana 
solutions exist in this case. 
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IV. MAJORANA SOLUTION FOR VORTEX IN THE SPIN-ORBIT COUPLING. 



For a planar system, the Rashba spin-orbit term a{cr x p) ■ z in the Hamiltonian we 
considered can also be written as acr ■ p. These 2 terms are simply related to each other by 
a az spin rotation and a more general Rashba-type spin-orbit term can be written as 

Hso = a[cosC(cr x p) ■ z + sm(cr ■ p]. (30) 

A recent proposaP^ has considered a defect in such a spin-orbit coupling where the angle of 
the spin-orbit ( varies in space to form a vortex (C(^) = 0)- The full BdG Hamiltonian for 
such a vortex can be written in Nambu space as 

HBdG = (-^V^ - /i)r, + Vz(Tz + ^ {a+ {a{r)e'',p^} + a_ {a{r)e-'' ,p+}) + At, (31) 

where the anti-commutation must be introduced in the spin-orbit term to preserve hermitic- 
ity. 

Substituting the circular-polar form for the derivatives we note that the Hamiltonian 
becomes 

1. 92 



HBdG = -Mdr + ^-dr " ^) + ^A^z + ^^^^ 

- '-{a^m. «W} + ^{e~*', e^'de}) + h.c} + At, 



(32) 



which in turn simplifies to a independent form 



HsdC 

1 



-a 



r 4^2^ r^i]Tz + Vza, 



ia{r) 



^{a+{{a{r)dr + a\r)] ^ 



{de + i}) + h.c} + At,. 



(33) 



Therefore, Jz = Lz commutes with the above Hamiltonian and the spinor form is 



raj \ 



(34) 



\ -V^,mA'r) I 

As before only the mj = channel can lead to a non-degenerate Majorana solution, and 
the BdG equation in this channel is given by 

1 



Hsdo = -{vid^ + -dr) + fJ-jrz + VzGz 
r 

- \{a+{{a{r)dr + a\r)} + ^) + h.c} + At, 
1 r 



(35) 



15 



^o(r) = 0. (36) 



(37) 



Since the above BdG equation is real, it can be reduced to a 2 x 2 matrix differential 
equation: 

^{^l + l^r) + V,-^l \A + a{r)dr + a'{r) + ^ 
-AA + a{r)dr + a'{r) + -^{^dl + ldr)-V,- ^^ 

As before, considering a step function vortex profile a{r) = for r < R and a{r) = a 
for r > R, one notices that the reduced BdG equation outside the spin-orbit vortex core 
resembles the reduced BdG equation in the same region for a regular vortex in the large r 
limit (Eq. ( [l5| with A(r) = A). Inside the vortex core a{r) = 0, and the BdG equations, 
as before, are analytically solvable via Bessel functions as below 

7]{-d^. - Idr) + V,-fi AA \ / u^Mzr) 

-AA r]{-d^, -ldr)-V,-fij [ u^Mzr) 

where 

-77^2 + V; - /i AA 
-AA -77^2 - 14 - 

As before, this leads to 2 solutions inside and 3 solutions outside the vortex core, with 5 
constraints at the interface. This leads to a single non-degenerate Majorana solution at the 
interface. 

Contrary to the result in Ref. [5S] for the asymptotic (r — )■ oo) behavior of the zero en- 
ergy wave function, il)Q{r) oc e~^^^/", we find that the asymptotic zero energy wave function 
behaves as il)Q{r) oc e~^'''/°' ^ for yU = and small A < V^. Therefore, according to our 
result, the decay length of the zero energy wave function diverges in the limit of vanish- 
ingly weak superconductivity (A — )■ 0) and the Majorana mode disappears by delocalizing 
over the entire system. This is in contrast to the result for the wave function in Ref. [HH] 
where the zero energy Majorana solution remains localized for arbitrarily small values of the 
superconducting gap. 




0. (38) 



V. MAJORANA SOLUTION ON THE SURFACE OF A TOPOLOGICAL INSU- 
LATOR 

Now we apply a similar argument to the vortex in proximity-induced s-wave supercon- 
ductivity on a TI surface 1^ which is obtained from our Rashba model by setting r/ = 0. The 
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equation for the allowed values of z in the superconductor for r > R are then 

/i^ - + {za ^ A)2 = (39) 
z = ±{A±t^/JF^VJ)/a (40) 

Therefore, in each of the ± channels, for small V^, there are a pair of complex conjugate 
eigenvalues on the same side of the imaginary axis. For the + channel both the eigenvalues 
are to the right of the imaginary axis and therefore are acceptable decaying solutions. Thus 
there are 2 linearly independent solutions (m-|-_i, w^^i) and (u-f-,2, M4„2) for r > R. From our 
previous discussion it is now clear that there is only one such solution {u^^3,ui^s) for r < R. 

Since the Hamiltonian is linear in the derivative, the boundary conditions only require us 
to match the wave-function (■u-)-(r), 'u^(r))^ at r = i? and not the derivative. The boundary 
conditions that the zero energy solution must satisfy at r = i? are given by, 

„ \fj««wy (41) 

n(R) ) \ ) U \ hAR) ) 

Together with the normalization condition for the zero energy wave function, the above 
equations provide 3 constraints for the three variables 01,02,03. This yields a unique zero 
energy Majorana wave-function for a vortex on a TI surface. 

VI. NUMERICAL CALCULATION OF THE VORTEX EXCITATION SPEC- 
TRUM IN THE SPIN-ORBIT COUPLED SEMICONDUCTOR. 

In previous sections we calculate and show the existence of a Majorana mode in a vortex 
at the interface of an s-wave superconductor and a spin-orbit coupled semiconductor. The 
most important information missing from these analytical calculations is the excitation gap 
above the zero-energy Majorana state, the so-called minigap. A proper calculation of this 
requires a numerical solution of the vortex problem which can be done by considering the 
system on a sphere with a vortex-antivortex pai j^^l^^ l as shown in the inset of Fig. ([3]^a)). 

The BdG Hamiltonian of this problem can be written as 

H = [r]p^p +^{{a X R) -p + p^ ■ {a X R)} - f,]T, + Vza ■ R + A(r)r, (42) 

where p = — ^[V — -R(i?- V)] is the non-Hermitean gradient operator restricted to the surface 
of the sphere and i? = ^. The above Hamiltonian takes a more familiar form in angular 
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coordinates as 

H = [^L^ -^L-a- /i]r^ + Vza ■ R + A{e) {cos (j)T^ + sm(f)Ty} 

= [^L" - ^{L.a, + \L+a. + \L_a+} - ^]r, + Vz{R.a, + + ]^R^a+] 

+ \^{R^r. + R_T^} (43) 

where R^ = cos 6* and R± = sin^e^**^. In these equations i?_ = R\_ and L_ = L^. The 
spectrum of excitations of this system is found by solving the eigenvalue problem 

H"^ = E^. (44) 

Similar to the vortex in the planar geometry, the BdG Hamiltonian has a combined 
spin-orbit-pseudospin rotational symmetry. This symmetry can be expressed compactly by 
noting that the Hamiltonian commutes with 

Jz = L, + ^{a,-T,), (45) 

where we have used the identity [R±, L^] = ±R±. The dependence of the eigenstates with 
mj = m can be written as e*^^'^\E'm(6') = e*(™"*^'^^"^^)/^)'^^m(6'). The independent part of 
the eigenstate \E'm(^) then satisfies a 1 dimensional BdG equation Hm^rniG) = i?m\E'm(^) 
where 

= Ul{<P)HU^{<P) (46) 

which can be explicitly checked to be independent, and Um{(p) = e*^'^"*^"'^"^"'-'/^-''^. To 
solve for the 9 dependent part of \E'm(6') it is necessary to convert Hm to a discrete matrix 
by expanding '^rn{9) = ^iQ,mPi^\cos6) where pI"^\cos6) are the associated Legendre 
polynomials which are the independent parts of the spherical Harmonics. In the associated 
Legendre polynomial basis the kinetic energy term has the simple diagonal form l{l + 1). 



Under the transformation in Eq. (46), the terms Rz,± in H transform into P|°'^^^(cos( 



in Hm- Therefore, its matrix elements in the associated Legendre polynomial basis can 
be calculated using the spherical harmonic addition theorem. A similar procedure can 
be used to calculate the matrix elements of the 6 dependent vortex. For a vortex, we 
take A(0) = Atanhi?sin^/^ where ^ is taken to be the length-scale of the vortex. From 
symmetry properties it is clear that A{6)/ sin 6 is an even polynomial in sin^ and can be 
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FIG. 3: (Color online) (a) Plots for the individual components of the 4-component wave-function 
^'(r) = (u-|-(r), ti^(r), t'j^(r), — w-f-(r))-^ for the zero energy Majorana state at the north pole. The 
components of ^ satisfy u„ = v* confirming the Majorana nature of the wave-function. We also 
show the semiconductor heterostructure on the surface of a sphere with a vortex and an antivortex 
(with reduced superconducting amplitudes at the vortex cores) situated at the north and the south 
poles. (b)Numerical results for the vortex mini-gap AE (solid line) and bulk gap (dashed line) 
plotted against the spin-orbit coupling strength a on the semiconductor. In these plots we have 
used A = 0.5, fi = 0.0 and tj = 1.0 in the units where Vz = 1- The spin-orbit coupling strength 
a = 0.3 in (a) and varies for (b). In these units, the size of the vortex core has been taken to be 
unity. 

written in terms of associated Legendre polynomials as 

^ = E^(2m)i^(Wi)(cos^) (47) 

where the associate Legendre polynomial can be written as P/(cos^) = — sin 6 P'l (cos 6) and 
Q = 2^(1^1) I-i Pi^{x)A{x)dx. As with the R operators, the angular momentum matrix ele- 
ments can be calculated from the above expansion by using the spherical harmonic addition 
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theorem. 

As in the analytic solution for the vortex, the angular momentum index m transforms 
from m — )■ — m under the particle-hole transformation S and we expect non-degenerate 



E = Majorana solutions of Eq. 44 only in the m = channel. This is confirmed by our 



numerical solution of Eq. |44] where we find that only in the topological phase Co < 0, are 
there a pair of states in the m = angular momentum channel whose eigenvalues approach 
exponentially with increasing radius R. The non-zero energy eigenvalue of the Majorana 
fermion is a result of the presence of 2 vortices in our calculation with a finite distance 
between them. The wave-function of the E = eigenvalue of the m = angular momentum 
channel, that is localized at the north-pole is plotted in Fig. (|3](a)). The components of the 
wave-function are seen to satisfy = v* confirming the Majorana character of these states. 
In the figure, the wave-functions of the Majorana modes is seen to decay and oscillate away 
from the North pole. The splitting between the Majorana modes into a pair of exponentially 
small oscillating eigevalues is a result of the overlap between the Majorana modes at the 
two pole^. 

Aside from the E = eigenvalue in the m = angular momentum channel, an isolated 
vortex confines a set of non-zero eigenvalues in other m ^ angular momentum channels. 
Of these, the eigenvalue with the smallest absolute value occurs in the m = 1 angular 
momentum channel and has an eigenvalue equal to the so-called mini-gap of the vortex. 
As mentioned before, the superconductivity in the non-Abelian superconducting phase is 
re-entrant with a bulk gap that is proportional to the spin-orbit coupling strength. As seen 
from Fig. (|3](b)) both the bulk and mini-gap are proportional to the spin-orbit coupling 
strength a. For spin-orbit coupling a ~ 1 and chemical potential /i = 0, both the mini- 
gap and the bulk gap are of order the induced pairing potential A. Therefore for the 
semiconductor structure where superconductivity is proximity induced, the mini-gap of a 
vortex can be tuned to be of order A if the chemical potential fi can be tuned to be less 
than order Vz, the Zeeman potential applied to the semiconductor. This is different from 
the conventional case of a regular (not proximity induced) where the chemical potential is 
of order Ef and the mini-gap is of order A'^/Ep which is much smaller than A. 

Thus the ability to independently control the chemical potential in the semiconductor 
heterostructure provides us with a powerful tool that can increase the mini-gap of the vor- 
tex in the semiconductor heterostructure shown in Fig. ([T]) by orders of magnitude from 
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the values in chiral p-wave superconductors to of the order of 1 K. This leads to the possi- 
blility of performing TQC with the Majoranas trapped in vortices in the heterostructure at 
temperatures which are as large as a fraction of a IK. 



VII. BULK TOPOLOGICAL QUANTUM PHASE TRANSITION. 

We found that the Majorana modes exist for a spin-orbit coupled semiconductor system 
only in the parameter regime Co = Kj-(A2 + /i2) > 0. This in turn was related to the parity 
of the roots in one half of the complex wave- vector plane of solutions outside the vortex core. 



As pointed out before in the context of Eq. (25), these roots are indeed properties of the 



reduced bulk superconducting Hamiltonian in the absence of a vortex. We also expect such 
a connection between the bulk properties and the existence of Majorana modes on general 
topological grounds 

To show explicitly the connection between the condition for the existence of Majorana 
modes (Co < 0) and the bulk properties, we note that even though the gap in the bulk 
superconducting state prevents the existence of propagating states at = 0, it allows 
evanescent states. Since the states at E = are particle-hole symmetric eigenstates \l/o of 



a real Hamiltonian, we can apply an argument analogous to Eq.(15) to obtain a bulk BdG 
equation in a A channel 



Considering an evanescent state of the form \l/o(a;,?/) = g-^C^'cose+j/sme)^^ leads to a con- 



straint on z which was previously written as Eq. (25). Therefore, the condition on Cq, which 
determines whether the phase supports a Majorana solution or not is precisely related to 
the parity of decaying evanescent modes in a given A channel in the bulk superconductor at 
E = 0. 

A change of the parity of the decaying evanescent modes requires an £^ = mode to 
become propagating, which can only exist if the bulk superconductor is gapless. Therefore, 
a change of the sign of Cq, which determines the topological nature of the phase, must 
be accompanied by a closing of the bulk spectrum. This is determined by the full BdG 
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Hamiltonian for a state with momentum fc(cos^, sin^) and can be written in the Nambu 
space as 

HBdG = ivk^ - + V,a, + ^(e-*V+ - e*V_)r, + Ar,. (49) 

The spectrum is obtained by considering Det^Hsdc — Ek) = which can be simphfied to 
the equation, 

El = + A^ + + a'k^ ± 2^V^/S.-^ + P{V^ + a2p) (50) 
where e = r^/c^ — /i. Setting /c = 0, it can be seen that 

El = {Vz ± v^A^T^)' (51) 

which vanishes as Cq becomes zero, as expected. Recent workP^ has shown that the quantity 
Co is the Pfaffian of the BdG Hamiltonian at = 0, (Co = Pf{HBdG{k = 0)ayTy)). The 
sign of Co, which determines whether the phase of the superconductor is non-Abehan or not 
has been shown to be relatecP^ to the parity of the first Chern number topological invariant 
describing time-reversal broken topological superconductors.'^^^ 

The phase diagram of the spin-orbit coupled semiconductor system can be understood 
from Fig. (|4]), which gives the variation of the quasiparticle gap versus the Zeeman splitting. 
One knows from topological stability of the Majorana fermion mode that, due to its non- 
degeneracy, the Majorana state is protected as long as the bulk gap does not close as 
one moves through the parameter space. In Fig. (|4]), the gap closes (at the wave vector 
k = 0) for the Zeeman splitting corresponding to = V^^ = A^ + yU^. The phase with 
Vz > Vzc supports the non- degenerate Majorana state while the phase with < Vzc does not. 
These two regions are separated by a gapless point in the parameter space, which signifies 
a topological quantum phase transition. The quantum phase transition is topological since 
the superconducting order on both sides is explicitly s-wave and the two phases differ only 
by the topological properties such as Majorana modes in defects and boundaries. 

VIII. COMPETITION BETWEEN SUPERCONDUCTIVITY AND ZEEMAN 
SPLITTING. 

The proposal to realize Majorana fermion modes in spin-orbit coupled semiconductor 
system involves the introduction of a large Zeeman potential. In general, a Zeeman splitting 
is known to compete with and eventually destroy superconductivity. To understand better 
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FIG. 4: Quasiparticle gap versus Zeeman coupling for various values of spin-orbit interaction a. 
The strength of the spin-orbit coupling in the inset is such that a = 0.3 corresponds to 0.1 eV- 
A. The gap vanishes at the critical value Vz = \/ + fj?. The spin-orbit coupling has negligible 
effect below this critical point and the superconducting gap is of a conventional s-wave type. Above 
the critical value and in the absence of spin-orbit coupling, the superconducting gap is destroyed 
by the Zeeman coupling. Spin-orbit coupling opens up a gap in this phase leading to re-entrant 
superconductivity with is topological. 

the competition between the Zeeman splitting and superconductivity in a spin-orbit coupled 
semiconductor, we first consider the case without spin-orbit coupling. This case is described 
by the BdG Hamiltonian 

HBdG = ivk^ - /i)r, + Vzcr, + Ar^. (52) 

The dispersion relation of this Hamiltonian is = ±14 ± + P. In this case, with 
Vz = 0, we obtain a conventional proximity induced s-wave superconductor with no Ma- 
jorana phase. As Vz increases above yZ/i^ + the quasiparticle gap of the system closes 
and one obtains a metal with a Fermi momentum kp given by rjkp = yU ± ^/V^ — A^. 
This is the well-known Chandrasekhar-Clogston limitP^ where strong Zeeman splitting sup- 
presses the superconducting quasiparticle gap. This suppression is due to the fact that, 
in the spin-polarized regime {\Vz\ > a small pairing potential cannot open a s-wave 
superconducting gap since the latter couples opposite spins. 

The BdG Hamiltonian ai kp is doubly degenerate and is given by 

HBdG = VW^rz + Vaz + At,. (53) 
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The degeneracy of the above Hamiltonian at the gapless point, which arises from the particle- 
hole symmetry, is lifted by the Rashba spin-orbit coupling term akpa^T^ in the semiconduc- 
tor. This yields a topological superconductor with a gap given by 

Eg ^ 2akF^. (54) 

Considering the g function of the Zeeman splitting (Fig. Q), it is clear that, 

for Zeeman splitting below the critical value Vz < a/A^ + /i^, the superconductivity is 
non-topological in nature. The topological superconducting phase that supports Majorana 
fermions is created by the application of a Zeeman splitting to suppress the conventional 
pairing potential. In this regime, the spin-orbit coupling can open up a gap resulting in a 
re-entrant superconducting phase. However, as is evident from the previous discussions, the 
re-entrant superconductivity is unconventional (topological) in the sense that it supports 
Majorana fermions. 

IX. TOPOLOGICALLY PROTECTED EDGE STATES AT INTERFACES. 

One of the signatures of a topological phase is the existence of gapless edge states 
which are inextricably linked to bulk topological properties such as Majorana modes in 
vortices.SESHSa The spin-orbit coupled semiconductor structure introduced in Sec. II can 
be shown to have gapless edge states using methods similar to the ones described in the 
previous sections. Furthermore it turns out that this approach to analyze the existence of 
Majorana edge modes does not impose additional requirements such as rotational invariance 
that were critical for the demonstration of a Majorana solution in a vortex. Therefore, this 
method can be used to examine the question of the existence of Majorana edge modes even 
in the heterostructures with more general forms of spin-orbit coupling proposed by AliceaP^ 
where the Zeeman splitting can be introduced by an in-plane magnetic field. 

A. BdG Hamiltonians for edges. 

Edges can be created in the semiconductor heterostructure by varying a parameter of 
the Hamiltonian such as /i, Vz or A perpendicular to the edge of a surface. Without loss 
of generality we can consider an edge that is perpendicular to the direction y. Because of 
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translational symmetry along the edge, the resulting edge BdG Hamiltonian has kx as a 
parameter. The momentum parameter transforms as k^ — > —k^ under the particle-hole 
transformation. Therefore, a non-degenerate Majorana mode can only exist for k^ = 0. 

Fixing kx = ^ reduces the two dimensional edge problem to a 1-dimensional BdG Hamil- 
tonian for a single band semiconductor with spin-orbit coupling (assumed to be linear in the 
momentum ky), which in general can be written as 

HBdG = {-Vdy - fJ^{y))rz + V^a- ■ B + ladyp ■ ar^ + A{y) cos (pT^ + A{y) sin (pTy (55) 

where the unit vector B is the direction of the effective Zeeman field and the unit vector p 
characterizes the spin-orbit coupling. Using the spin rotation transformations on Hbcig, we 
can choose p — y without loss of generality. This yields the Hamiltonian 

HsdG = {-r]dl - iJi{y))Tz + V^ct ■ B + a{iay)Tzdy + A{y) cos (pr^ + A{y) sin (pTy (56) 

which is invariant under spin-rotations about the |/-axis. Therefore without loss of generality 
we can reduce the above Hamiltonian to 

HBdG = {-Vd^ - fJ'iy))Tz + Vz{cos + siu vGy) a{iay)Tzdy A{y) cos (t)Tx A{y) sin (t)Ty. 

(57) 

Non-degenerate Majorana spinor solutions are of the form ^ = {u,tayU*) and are com- 
pletely determined by the 2- spinor u. This fact was used to obtain the Majorana solutions 
for vortices to reduce the BdG equation from a 4 x 4 system of equations to a 2 x 2 system 
of equation. However, this reduction procedure required the BdG Hamiltonian to be real 
which is not the case for general forms of spin-orbit coupling and Zeeman splitting. The 
BdG equation for the zero energy mode HBdc^ — may be reduced to an equation for u as 

[{-dl - n{y)) + Vz{cos + sin voy) + a{iay)d^ u + A{y)e"^{%ay)u* = 0. (58) 

This equation is not real but may be reduced to a system of real equations by writing 
u = uji^ lUj and taking the real and imaginary parts of the resulting equation giving a pair 
of equations of the form 

\{-dl - H.{y)) + Vzco^voz + a{zay)dy + A{y) cos (f){zay)] ur - [A(y)sin0 - Vzsmi']{zay)ui 

(59) 

[{-dy - fJ>{y)) + Vz cos ua^ + a{iay)dy - A{y) cos (t>{icry)\ uj [A{y) sin - Vz sin v\{iay)uR 

(60) 
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The above pair of equations is similar to the pair of equations obtained for the two A = ±1 
channels except that earlier the two channels were decoupled. For A{y) independent of y, 
the 2 channels can also be decoupled by choosing such that A{y) sin (f) = Vz sin u. In what 
follows, we will make this choice and also replace Vzcosu — > Vz and Acos0 — )■ A. This 
results in a reduced BdG equation for the E = reduced spinor \l'o(y) 

y -AA(!,) - ad, -Tidl -V,- My) j 

where A = ±1. 

An edge in a two dimensional system of the type considered above is defined by requiring 
some parameter of the Hamiltonian to vary across the edge situated at ?/ = 0. We take this 
parameter to be constant for ?/ < and ?/ > 0. In this case, our previous approach can be 
applied in a way even simpler than the application to the vortex problem, since the solutions 
on both sides of the interface at ?/ = can be approximated as a sum \l/o(y) = Xln '^n^~^"^'^n 
where, as in the interior of the vortex (but far from the vortex core), (Eq|24)), 

-rizl + Vz- AA - Zna 
-AA + Zna --qzl - - fi 

Similar to the vortex case, in the topological phase Co = (A^ + /i^ — V|) < 0, there are 
3 values of Zn such that Re{zn) < 0, while in the non-topological phase Co > 0, there are 
only 2 solutions in a given A channel. The coefficients a„ in the solution are determined by 
matching the boundary conditions on \I^o(?/) at y = 0. The coefficient Cq, written in terms 
of the original parameters of the wire, reduces to 

Co = A^ cos^ + /i^ - cos^ z/ = A^ + /i^ - (63) 

and is not affected by the (p and u parameters that were introduced to make the BdG 
Hamiltonian real. The procedure of reducing the BdG Hamiltonian to a real Hamiltonian 
only introduces the additional constraint |Vzsinz/| < A. 

B. Chiral edge states. 

Based on analogy with FQHE and chiral p-wave superconductors, one expects a chiral 
gap-less state confined to the edge of the semiconductor heterostructure. An edge can be 
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created in such structures by raising the chemical potential /i towards the edge such that 
electrons stay confined inside the system. Therefore, an edge of a system confined to y < 
is defined by fi{y) = ior y < and fj, (y) = oo > IV2I for |/ > 0. In these structures we 

assume that A(y) = A is independent of y. 

The BdG equation now reduces to a 2 x 2 system of equations, 

[-fjdl -11 + Vza^ - lauydx + i\ay^]ilj{x) = 0, (64) 

where 

i^ix) = I '''^''l I . (65) 
In order to solve a semi-infinite system we make a plane-wave trial solution ansatz 



iP{x) = e"" ' (66) 




where z must now satisfy 

—rjz^ + Vz — n za \ ( 



= 0. (67) 



az -rjz"^ -V^- 11 j \ui 

As in the case of the vortex, for > (A^ -|- /x^) there are 3 solutions on the right 
half of the complex z plane and 1 solution on the right half for A = —1. The situation is 
opposite for A = 1. Solutions with Im{z) > are physical on the left-edge of the system 
while Im{z) < is physical on the right edge of the system. Thus for A = —1 there are 
3 physical solutions on the left-edge of the system which is the exact number needed to 
make the 2-component spinor vanish at the left-edge. Consequently, there is a localized 
zero mode at the left edge of the system in the A = — 1 channel. Similarly there is a 
locahzed zero mode at the right edge of the system in the A = channel. Finally, for 
< (A2 + there are no zero energy solutions at either edge since there are only 
2 physical solutions at each edge which is insufficient to match the boundary conditions. 
Since the wave-function is confined to the edge, we expect the boundary conditions ^o,t(0) — 
^0,4.(0) = 0, which together with normalization lead to 3 constraints. As mentioned before, 
in the topological phase we obtain 3 coefficients corresponding to the 3 normalizable 
solutions in the interior. Therefore, there is a unique zero energy state resulting from a 
matching of the boundary conditions. This state is a Majorana mode for the end-point of 
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the nanowire in the topological phase(Co < 0), which disappears when we tune the wire 
through a phase transition to Cq > 0. The Majorana modes at the edges discussed above 
only occured at = 0. The complete spectrum of the edge is obtained by considering the 
BdG Hamiltonian at small ^ using the k ■ p perturbation theory. To lowest order in 
kx, the Nambu Spinor wave-function can be approximated as \1/A;^(a;, ?/) ~ e*'^"'^\l/o(?/) with 
an energy Ek^ = vk^ where v = (\l/o|cryr^|\l/o)- 

A similar chiral Majorana wire is obtained by considering an edge between the topological 
phase Cq < and Cq > where fi{y) is constant but Vz{y) = for ?/ > 0. In that case there 
are 5 constraints to match as in the vortex case, and there are 5 coefficients, 3 arising from 
the topological phase at y < and 2 from the non-topological phase y > 0. 



C. Non-chiral Majorana edge states. 

Now we consider the junction of a pair of topological superconducting islands with phases 
(j) and (f)' = T[ — (f) which is a geometry that is of particular interest to TQC architechtured^^'^. 
For such a choice of phases, the effective pairing potential in A cos is a step function 
given by A(x) = A cos for x < and A(x) = — A cos cj) ioi x > 0. As before we then replace 
Acos0 — i- A. Focusing on the k^ = Q particle-hole symmetric momentum for the edge 



, solutions for X < and a; > can be expanded in terms of spinor functions given in Eq. 66 
which is written in terms of eigenvalues z and eigenvectors satisfying Eq. [67| Normalizable 
solutions must now be composed of superpositions of exponentials with Re{z) < for ?/ < 
and Re{z) > for x > 0. We note at this point that the equations for x > and x < 
differ by a change in sign of A(x) across the interface which corresponds to a change in sign 
of z. Thus, as before for A = —1, in the topological phase we have 3 values of z such that 
Re{z) > for x > and 3 values such that Re{z) < for x < 0. Following the boundary 
condition matching argument of the last section, for the vr junction there are 6 states at 
X = to compose wave-functions at x = which need to satisfy 5 constraints. Therefore, 
generically there will be a pair of zero energy modes satisfying these equations. 

It might appear that, unlike the case for the chiral edge states, the pair of Majorana states 
cannot be topologically protected. In the case of the TI/SC interface, |37j the existence of 
such a pair of non-chiral Majorana edge modes at a phase difference vr was a consequence of 
time-reversal symmetry which is broken here. In our calculation we find that because of this 
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time-reversal breaking, the pair of degenerate zero modes may occur at a phase difference 
of TT — 2(j). In fact, by considering the evolution of the Andreev bound state spectrum in 
the junction as a function of phase difference, ^ it is possible to show that even though the 
Majorana nature of the pair of non-chiral Majorana modes is not topologically protected, it 
is not possible to eliminate the zero-modes all together. The zero-crossing of the non-chiral 
Majorana modes may only be shifted to different values of phase by time-reversal breaking 
perturbations. This is a result of the fact that the 2 particle-hole symmetry related branches 
of the Andreev bound state spectrum differ by Fermion parity. Therefore, even though an 
infinitesimal perturbation can lift the degeneracy of the 2 zero modes at some value of phase 
difference 50, it can only do so by shifting the crossing to a neighboring value of 5(f). Similar 
to the chiral edge modes one can use the k ■ p perturbation theory to construct a pair of 
linearly dispersing modes from the pair of zero energy states at ky = 0. 

X. MODEL CALCULATIONS OF PROXIMITY EFFECTS IN SUPERCONDUC- 
TOR - SEMICONDUCTOR - MAGNETIC INSULATOR HETEROSTRUCTURES 

In this section we study, starting from a microscopic tight binding model, the excita- 
tion spectrum of a semiconductor thin film sandwiched between an s-wave superconductor 
(SC) and a ferromagnetic insulator (MI). We determine the dependence of the effective SC 
gap and Zeeman splitting induced by proximity effects on the parameters that characterize 
the heterostructure model. We also calculate the dynamical contributions to a low-energy 
effective theory of the proximity effect and identify parameter regimes suitable for the exper- 
imental implementation of the semiconductor-based proposal of a platform for topological 
quantum computation. 

To study the proximity effect in the SC - semiconductor - MI heterostructure, we consider 
the minimal microscopic model defined by the Hamiltonian 

H,,t = Ho + Hsc + Hmi + Hi, + H-,^^ . (68) 

The Hfj term describes the semiconductor thin film, 

i,j,cr i,cr,cr' 

- cl^5yai^^x)ayCi^' + h.c] + ^ V{zi)clci^, (69) 
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where the first contribution describes hopping on a cubic lattice, while the second represents 
a lattice model of the Rashba spin-orbit interaction. The hopping matrix elements are non 
vanishing for nearest neighbors, tjj = —to? and next nearest neighbors, tjj = —ti, and we 
also include an on-site contribution ta = eo that shifts the bottom of the semiconductor 
spectrum to zero energy. The parameter a represents the Rashba coupling constant and 
6x(y) are nearest neighbor displacements in the xy-plane of a cubic lattice with lattice pa- 
rameter a. The system is assumed to be infinite in the x- and y-directions and contains N 
planes perpendicular to the z-direction. The quantities f^i^y) are Pauli matrices and c|^, Ci^ 
are electron creation and annihilation operators, respectively. The last term represents an 
external bias potential that modify the on-site energies along the z-direction. 

For the superconductor we use a simple mean-field model defined by the Hamiltonian 

Hsc = -ts Yl blbja + ^ blh^ (70) 



where ts represents the nearest neighbor hopping on a cubic lattice and A is the mean-field 
s-wave SC order parameter. The SC and the semiconductor thin film have a planar interface 
perpendicular to the z-direction. For clarity, the electron creation and annihilation operators 
inside the SC were denoted bl^ and bi„, respectively. 



The third term in Eq. (68 ) represents the ferromagnetic insulator which, again, is modeled 



at the mean field level by the Hamiltonian 

ataia. (71) 



HmI — — E] tmaO'la'^jiT + E/ f^""^ ~ ^ 2^ ] 



The nearest neighbor hopping is spin dependent and we consider the case = — ^mt > 0- 
The spin-dependent on-site energy was divided into a contribution = Cmo + Qtma that 
places the top of the valence band and the bottom of the conduction band at the same 
energy e^o and a term proportional to the insulating gap F. The last two terms in Eq. 



(68) describe the coupling at the two interfaces, semiconductor-SC and semiconductor-MI, 



respectively: 

^^s = - E i5(&.W + h.c.), (72) 
Hh. = - E iM(aLc,. + h.c.). (73) 
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The parameters ts and Im characterize the transparencies of the two interfaces and provide 
the energy scales for the couphng between the semiconductor thin film and the SC and MI, 
respectively. These energy scales are crucial for determining the strength of the proximity 
effects. 



We diagonalize Eq. (68) numerically for a system with periodic boundary conditions 
along the x- and y-directions and a finite size along the z-direction. The semiconductor, 
SC and MI contain N, Ng and Nm planes, respectively. For the semiconductor thin film 
we consider 5 < < 20, while A^^ and Nm are typically of the order of several hundreds. 
For large systems, the proximity effects are independent on A^^ and and we explicitly 
checked that the values used in the calculations are within that regime. The hopping matrix 
elements for the semiconductor are to = ^-^TeV and ti = 0.39eV^, which, for a cubic lattice 
with a = 5.5A generate in the low wavelength limit a small effective mass m* = 0.04me 
characteristic of semiconductors such as InAs. The value of the Rashba coupling constant is 
a = 18 meV, corresponding to a strong spin-orbit coupling of the order 100 meV A. For the 
SC, the hopping parameter is ts = 0.35e\^ and A = 1.5meV, while the magnetic insulator is 
described by t^i = —tm-\ = 0.36eV^ and a gap value F = 200meV. The coupling parameters 
ts and tM are varied and typically range within several hundreds meV. 



A. The magnetic proximity effect 

First, we turn on the coupling at the interface between the semiconductor and the MI. 
Fig. [5] shows the spectrum of a slab containing a semiconductor thin film with A^ = 10 
atomic layers in contact with a MI. The coupling at the interface is tm = 250meV. The 
semiconductor spectrum is characterized by 10 strongly dispersive bands represented by red 
curves in the upper panel of Fig. [5j For an isolated semiconductor, these bands are weakly 
split by the spin-orbit interaction and are double degenerate at k = 0. However, as a result 
of the magnetic proximity effect, an effective Zeeman splitting removes this degeneracy. This 
is shown in the lower panel of Fig. [5] for the lowest energy mode that has the minimum 
inside the insulating gap. The induced Zeeman splitting is entirely a result of the exchange 
interaction between the MI and the SM layers and is not related to the magnetic field 
produced by the ferromagnetic MI. In fact, effective Zeeman splittings have been known to 
induced by anti-ferromagnetic insulators as well through a phenomenon commonly known 
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FIG. 5: (Color online) Band structure for a semiconductor-ferromagnetic insulator heterostructure 



described by equations (68), (69), (71) and (73) with = 10 and tM = 250meV. The red points 
represent semiconductor states, while the magenta bands are ferromagnetic insulator states. The 
lower panel shows the low-energy behavior around k = 0. Notice the spitting of the semiconductor 
band at k = due to the effective Zeeman term induced by ferromagnetic proximity effect. 

as exchange bia^. 

A natural question of practical importance concerns the dependence of the induced Zee- 
man splitting on the parameters of the model. To get a better intuition of the physics 
behind the proximity effect, we start with an analysis of the structure of the wave functions 
of the relevant states. Consider a semiconductor state with an energy within the insulat- 
ing gap r. When the semiconductor and the MI are coupled, the wave function describing 
this state will partly penetrate into the magnetic insulator, where it decays exponentially. 
These components of the wave functions inside the MI will acquire ferromagnetic correla- 
tions equivalent to an effective Zeeman field. Hence, the strength of the proximity effect is 
related to the fraction of the wave function that penetrates the MI. This qualitative picture 
of the proximity effect is illustrated in Fig. |6} In general, the fraction of the wave function 
that penetrates the insulator depends on three types of factors: i) the properties of the 
quasi-two dimensional parent system, in this case the semiconductor film, ii) the properties 
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FIG. 6: (Color online) Dependence of the amplitude of low-energy states on the distance from 
the semiconductor - MI interface for three different sets of parameters. In order of increasing 
penetration depth into the MI: lowest semiconductor band (n = 1) and tm = 0.25eV (blue), second 
semiconductor band (n = 2) and tm = 0.25eV (black), and n = 2, tm = 0.5eV (red). For a given 
set of parameters, a certain fraction of each low-energy wave function penetrates the insulator 
leading to ferromagnetic correlations in the semiconductor thin film,, i.e., to an effective Zeeman 
field induced by proximity effect. 

of the insulating host system, and iii) the coupling strength at the interface. 

The first and last types of factors are qualitatively illustrated in Fig. |6j For example, by 
controlling the value of the on-site energy eo, i.e., chemical potential of the semiconductor, 
one can bring the minima of different semiconductor bands within the insulating gap. The 
profiles of the wave functions in the transverse direction (i.e., perpendicular to the film) 
depend on which band is inside the gap, as shown in Fig. [6] for the first two bands, n = 1 
and n = 2. Furthermore, the fraction of the wave function that penetrates the insulator 
depends on its amplitude at the interface. This amplitude is significantly different for n = 1 
and n = 2, hence one expects significantly different strengths of the proximity effect. Of 
course, the amplitude at the interface can be also modified by changing the size of the system 
(N), or by applying a bias potential V{z) that can tilt the spectral weight toward or away 
from the interface. In addition to the amplitude at the interface, which is determined by 
the parameters of the semiconductor, the fraction of the wave function that penetrates the 
Ml also depends on the slope at the interface. In essence, this slope is controlled by the 
coupling parameter im, as illustrated in Fig. [g] Finally, the parameters of the exponential 
decay inside the insulator, as well as the strength of the ferromagnetic correlations depend 
on the properties of the insulating host. 
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FIG. 7: (Color online) Dependence of the normalized effective Zeeman splitting on the semiconduc- 
tor thickness for different sets of parameters (dots). The continuous lines represent the normalized 
wave function amplitude at the interface (the Nth layer) |'0„(A^)p/|^/^„(lO)|^. The reference gap 
is Tq = 2.955 meV. The wave function amplitude in the absence of an applied bias, V{z) = 0, is 



determined using Eq. (76), while for a linear bias V{z) = Vqz/w with Vq = ±1 eV the amplitude 



is determined numerically. 

Next we proceed with a quantitative analysis of the induced Zeeman splitting and derive 
an effective low-energy theory for the ferromagnetic proximity effect. We consider the sit- 
uation when the states of the nth semiconductor band with wave vectors near k = have 
energies inside the insulating gap. The states ipnkuij^i) corresponding to vanishing spin-orbit 
coupling, a = 0, form a convenient basis for the low-energy Hilbert subspace of interest. 
Projecting the Hamiltonian Hq onto this subspace we obtain 



a{ky — ikr, 



ffoeff(k)= I p h (74) 

a{ky + iK) 2^ 

where m* is an effective mass with a value slightly different from the effective mass corre- 
sponding to the 3D semiconductor model. The difference stems from the quasi 2D geometry 
of the system. The insulating degrees of freedom can be integrated out and replaced by 
an interface self-energy. When projected onto the low-energy subspace, this contribution 
becomes 

S..,(k,a;) = -tl|^„k(^o)|'Gi™Hk,^; + 5.), (75) 
where G^J^]{k,u; zo + 6z) is the Green function of the magnetic insulator and |'?/'k(-Zo)P the 



amplitude of the semiconductor wave function, both at the interface. From Eq. (75) one 
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FIG. 8: (Color online) Dependence of the normalized Zeeman splitting on the interface transparency 
for a system with = 12 and different model parameters (dots). The straight lines are guides 
for the eye. The reference interface coupling is tmo = 250 meV. The rather small deviations from 
a linear dependence indicate that dynamical effects are negligible, i.e., neglecting the frequency 



dependence in Eq. ( 75 ) is a good approximation. 



immediately notice that, neglecting dynamical effects, i.e., setting u = 0, any induced 
Zeeman splitting has to be proportional with the amplitude of the wave function at the 
interface times the square of the interface transparency. For a vanishing bias potential, 
V{z) = 0, the wave function amplitude on the j atomic layer is proportional to sin^[7mj/ (A^+ 
1)] and the amplitude at the interface becomes 

, 2sin2(^) 

\A.uzo)\'= (76) 

Note that, for clarity, we dropped the spin label, as the amplitude is spin independent. Also 
note that for large N the amplitude becomes l^/'nkl^o)!^ ~ 2n^7r^/(A^+ 1)^. The dependence 
of the induced Zeeman splitting on the wave function amplitude is shown in Fig. [7j As the 
film thickness w = Na is increased, the value of the of the wave function amplitude at the 



interface drops rapidly, as expected from Eq. (76). Due to the proximity effect, an effective 



Zeeman splitting creates a gap in the nth semiconductor band at k = 0, 

f = E„2(0)-E„2(0), (77) 

where £'„j(k) are the energies of the nth semiconductor band within the magnetic insulator 
gap (see lower panel of Fig. [s]). The proportionality between the induced Zeeman gap f 
(dots in Fig. [T]) and the amplitude of the wave function (lines in Fig. [T]) reveal the absence 
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FIG. 9: (Color online) Dependence of the effective Zeeman splitting on the ferromagnetic insulator 
gap for two sets of semiconductor parameters (dots) . The continuous lines are calculated using the 



effective low-energy theory described by Eq. (81). Notice the remarkable agreement between the 



low-energy theory and the exact results and the difference in energy scales between the insulating 
gap and the induced gap. 



of significant dynamical effects Eq. (75). This conclusion is further supported by the linear 
dependence f on the square of the interface transparency, t^, shown in Fig. [sj We empha- 
size that the effective coupling constant that determines the strength of the ferromagnetic 
proximity effect, gm = 2t^|'?/'n(-2o)PMm, where is a characteristic bandwidth for the 
magnetic insulator, can be tuned by: a) modifying the semiconductor film thickness (see 
Fig. [T]), b) applying a bias voltage (see Figs. [7]and[8]), and c) changing the semiconductor 
- MI coupling (Fig. 

So far we have discussed the dependence of the proximity effect on the properties of the 
quasi-two dimensional parent system and on the coupling strength at the interface. Next 
we briefly investigate the role of the insulating host system. We note that a full analysis 
would require a treatment of the magnetic insulator beyond the simple mean-field picture 
used here, but this would not change the results obtained so far. Within our mean-field 



approximation, the MI Green function from Eq. ( 75 ) can be easily evaluated 



G;,'^)(k,a;;^o + 5.) = 5..'5^^Twi;v— -Ix.k.(^0 + 5.)P, (78) 

where E^f^{k) are the energies of the insulator bands and Xi'ktri^o + the values of the 
corresponding eigenstates at the interface. It is convenient to express the Green function in 
terms of the partial density of states Pko-(i^) = "^i^^i^ " E^^{k))\xuka{zo + 5^)P- Within 
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FIG. 10: (Color online) BdG spectrum of the full Hamiltonian (68) for a semiconductor thin film 
with = 10 sandwiched between a MI (F = 0.2eV, tm = 250meV, F = 2.95meV) and an s-wave 
SC (A = l.SmeV, ts = 130meV). The red points represent states that reside (mostly) within the 
semiconductor, magenta designates the MI bands, while the SC states are blue. Due to the effective 
Zeeman filed, the semiconductor bands are split and only one crosses the chemical potential (see 
upper panel). As a result of the superconductor proximity effect, a small gap opens at the crossing 
points (see lower panel). 



our model, this partial density of states becomes 



oo-E^ik) 



(79) 



with Ag- = 2|tmo-| being half of the bandwidth in the insulating phase and -E'o.(k) the energy 
values at the middle of the valence and conduction bands for a given wave vector parallel to 
the interface. At k = we have -E'o-(O) = — yU — cr(r/2 + Aq-). Note that pko-(i^) vanishes for 



values of tu outside the bandwidth. Using Eq. (79), we obtain for the Green function the 
expression 



2 



sign(cu - -Ect)^ 



ou-E^ 



- 1 



(80) 



A. ° ' V V A. 

Note that the imaginary part of the Green function vanishes for values of u within the 
insulating gap. Also, because the energies of interest are much smaller than the insulator 
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FIG. 11: (Color online) Dependence of the low-energy spectrum on the coupling between the 
semiconductor and the SC: ts = (red), tg = 130meV (black), tg = 190 (yellow). A proximity 
effect - induced gap opens at finite k. As tg is increased, the minimal gap shifts to lower wave 
vectors and first increases, then decreases and eventually vanishes at a critical coupling before 



opening again (see also Fig. 12) 



bandwidth, a; ^ Aq-, we can neglect the frequency dependence in Eq. (80). Finally, within 



the static approximation, we obtain for the induced Zeeman splitting the expression 



AF^ - AFt, 



(81) 



AF. 



-fi - (t(| + A^) 



\ 



-fl - 0-(| + Aa) 



where the effective coupling constant is gm = 2t^|^/'„(2;o)P/Ao-. To test the accuracy of 
this effective low-energy theory, we compare the values of the induced Zeeman splitting 



predicted by Eq. (81) with the numerical calculations. The results for various values of the 



insulating gap shown excellent agreement (see Fig. [9]). Finally, we note that away from 
k = the dispersion of the low-energy bands £'„j(k) can be obtained by adding the self- 



energy contribution (75) to the effective theory described by Eq. (74). Within the static 



approximation we have 



i^eff(k) 



2m* 2 0!{ky ikx) 



(82) 



with F given by Eq. (81 ). Explicit calculations for various sets of parameters show that this 



low-energy theory represents an excellent approximation for all k values of interest. 
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B. Superconducting proximity effect 



Next, we turn our attention to the effects induced by the proximity of an s-wave supercon- 
ductor on the semiconductor-MI heterostructure. The parameters of the semiconductor-MI 
interface are fixed, with tm = 250 meV. We consider a semiconductor thin film with = 10 
and create a new interface at the free surface of the semiconductor by couphng it to a SC 



with an s-wave gap A, i.e., we add the terms given by equations (70) and (72) to the total 



Hamiltonian (68). The corresponding spectrum is shown in Fig. 10 The semiconductor 
band is split due to the ferromagnetic proximity effect and the chemical potential is tuned 
so that it crosses only the lower energy mode. When the coupling to the SC is turned on, a 



small gap opens at low energies due to the SC proximity effect (see Fig. 10 lower panel). 
Before a quantitative analysis, let us illustrate qualitatively the behavior of the proximity- 



induced SC gap. Fig. 11 shows the low-energy spectrum for three different values of the 
coupling constant tg. For tg = (red line) the BdG spectrum is gapless. A non-vanishing 
interface coupling opens a small gap at the crossing points. The value of this finite wave 
vector gap increases with tg, but the gap at k = decreases, (black line). Eventually, at a 
critical value tgc the gap vanishes at k = 0, before opening again for larger couplings (yellow 
curve). This closing of the induced gap signals the presence of a quantum phase transition!^. 

In order to understand this behavior, it is useful to develop an effective low-energy theory 
for the SC proximity effect. As shown previously,^^^! dynamical corrections are crucial in 
capturing the low-energy physics in this case, in contrast to the ferromagnetic proximity 
effect. Using the results obtained in Refs. [I2|l3] . the Green function describing the low- 
energy physics of the heterostructure can be written as 

- = (a + \ka+ + A^a_)r, + ^a, + , - to ( 1 + ^ (83) 



2 VA2 - w2 ^ ^ 



where = k'^/2'm* — /i, F is given by Eq. (81 ), = a{ky — ikx), and the effective coupling 
is Qs = 2i'j.\ip{zs)\'^ / Ag. The wave function amplitude is evaluated at the semiconductor - 
SC interface and A^ = 2ts is half of the SC bandwidth. The low-energy spectrum can be 
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FIG. 12: (Color online) Upper panel: Comparison between the solution of the effective low-energy 



theory given by Eq. ( |84[ ) (lines) and the numerical solution of Hamiltonian (68) (dots) for two 
values of the chemical potential. Only positive energies are shown. Lower panel: Dependence of 
the induced minimum gap on the interface transparency. The vanishing of the gap at a critical 
value tscip-) reflects a quantum phase transition between a topologically non-trivial SC (at small 
values of is) and a trivial s-wave SC (at large couplings) [reference]. 



obtained by solving the corresponding BdG equation, Det(G' ^) = 0. Explicitly, we have 

2 f2 „2a2 



1 + 



9s 



VA2 - 



_Lf2 I |x |2 I 9s^ 



A 



+ 



4 A2 -a;2- 



(84) 



A comparison between the solution of Eq. (84) and the numerical calculations is shown in 



the upper panel of Fig. 12 The good agreement between the two calculations indicates that 
all the relevant ingredients have been incorporated into the effective low-energy theory. The 
dependence of the minimum gap on the interface coupling is shown in the lower panel of 
Fig. 12 Note that the critical value of ts at which the gap vanishes can be obtained from 
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FIG. 13: (Color online) Details of the low-energy spectrum for the case when the n = 2 semicon- 
ductor band has the minimum in the vicinity of the chemical potential. The low-energy physics in 



the vicinity of k = (upper panel) is described by the effective theory given by Eq. (84) with a 
coupling constant gs that includes the amplitude of the n = 2 state at the interface. The n = 1 
semiconductor bands cross the chemical potential at larger values of k (lower panel) and are gaped 
due SC proximity effect. Note that the gap at large wave vectors does not vanish at a critical 
coupling. 



Eq. (84) by setting a; = 0. Explicitly we have 




The last question that we address in this section concerns the situation when a higher 
energy semiconductor band, n > 1, has the minimum in the vicinity of the chemical potential. 
For concreteness, we consider the case n = 2. The physics in the vicinity of k = is similar 
with the case studied above. In addition, the n = 1 bands cross the chemical potential 
at some large value of k. Nonetheless, assuming that the partial density of states of the 
superconducting metal piiai^^) does not vary significantly with the wave vector, i.e., the 
effective mass of the metal is much grater than the effective mass of the semiconductor, the 
n = 1 bands will be gaped and the induced gap is typically grater than the minimum gap 
near k = 0. Hence, the general conclusions of our analysis of the n = 1 case hold for n > 1. 
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To illustrate this point, we show in Fig. 13 the relevant details of the low-energy spectrum 
for the case n = 2. 



XI. MAJORANA FERMION MODES IN ONE DIMENSIONAL NANOWIRE. 

In the previous sections, we discussed how Majorana states may appear at vortices and 
edges of various two dimensional spin-orbit coupled semiconductor heterostructures. In this 
section, we show that Majorana fermions can also be realized in the much simpler one- 



dimensional nanowire set up (Fig. 14). In this set-up we propose to study a semiconducting 
nanowire with a sizeable spin-orbit coupling ( e.g., InAs) placed on an s-wave superconductor 
{e.g.., Al). An in-plane magnetic field is used to create a Zeeman splitting in the nanowire 
which gives rise to the band-structure shown in Fig. [I4|(b). The direction of the magnetic 
field (parallel or perpendicular to the length of the wire) that is required to open a gap in 
the nanowire depends on the exact direction of the spin-orbit coupling in the wire. The 
chemical potential in the wire /i is assumed to be controlled by external gate voltages. 



We argue below that, for the Zeeman splitting satisfying > A/yU^ + A^, a single non- 
degenerate zero energy state exists at the end of the wire as the only low-energy bound state. 
The second quantized operator for this state is again a Majorana fermion operator. This 
mode can be detected as a zero-bias conductance peak in the STM tunneling spectrum. On 
reducing the Zeeman splitting so that Vz < a//!^ + A^, the Majorana mode, and hence the 
zero-bias conductance peak, should disappear from the tunneling spectrum. 

The BdG Hamiltonian for a one-dimensional single band semiconductor with spin-orbit 
coupling (which is linear in the momentum) in proximity to an s-wave superconductor (Fig. 



14), can be written as 



Hsdc = {-Vdy - iJ^{y))Tz + V^CT ■ B + ladyp ■ (tt^ + /\ cos ^r^, + A sin (pTy. (86) 

Here the unit vector B gives the direction of the effective Zeeman field and the unit vector p 
characterizes the spin-orbit coupling. From inspection it is now clear that this Hamiltonian 



is formally identical to the Hamiltonian for the edge states, Eq. 56, at = 0. Therefore, 
from the solution of the chiral edge state at = for the semiconductor thin film (Sec. X 
B), we conclude that in the topological phase of the wire (Co < 0), the ends of the nanowire 
support localized zero-energy Majorana states. 
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FIG. 14: (a): Geometry to detect zero-energy Majorana fermions using STM spectroscopy on a 
semiconducting nanowire. The Zeeman splitting is induced by a parallel magnetic field, while the 
chemical potential is controlled by an external gate (not shown). The Majorana fermion mode 
localized at the end of the nanowire gives rise to a zero-bias peak in the STM tunneling spectrum 
from the end. The tunneling spectrum from the bulk of the nanowire is gapped, (b): Band- 
structure of the nanowire in the topological superconducting phase. 

A. STM detection of Majorana end modes in the nanowire. 

In the previous section, we showed that the end of a nanowire in the topological phase is 
characterized by a Majorana mode. The Majorana modes at the ends of a one dimensional 
p-wave superconductor have been shown to lead to distinct signatures in the STM tunneling 
spectrumJ ^^ I ^^ I In what follows, we analytically and numerically calculate the STM conduc- 
tance spectrum from one end of our semiconducting nanowire. We find that, for the Zeeman 
coupling satisfying {Vz > a/AM-T?), this spectrum has a zero-bias conductance peak. The 
zero-bias peak disappears as the Zeeman splitting is reduced to satisfy (Vz < a/A^ + /i^). 

The tunneling spectrum of a superconducting system, similar to the tunneling conduc- 
tance of normal systems, can be calculated using the Keldysh formalism of non-equilibrium 
Green function^ where the superconducting system is coupled to a tip which is initially 
in thermal equilibrium at a chemical potential fit, by an adiabatically increasing tunneling 
amplitude perturbation V. To calculate the current, we consider an STM tip state ipKx') 
where x' is restricted to the tip. The tunneling Hamiltonian between the STM tip and the 
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superconducting wire can be written as 

Htunnel = I dxdx' J^l^ M^U^')Mx) + V* {xx')^l{x)Mx')] ■ (87) 

The tunneling current can be related to the Keldysh Green function of the combined system 



a^ 



I = f dxdx'J2[^*i^^')Gif{xx') - V{xx')Gi''J{x'x)]. (88) 

The Keldysh Green function G^^^ can be evaluated using the Dyson equation 

Gi^Axx') = -GfXxxr)V{x,x\)gZ\x',,x') - Gf:{xx{)V{x,x\)g^J^\x\,x') (89) 
GW(x'x) = -gi^\x'x\)V*{x,x\)GZ\x,,x) - gif{x'x[)V*{x,x[)Gii\^i,^) (90) 

where g = {Hq—uj) are the unperturbed Green functions and Hq is the unperturbed Hamilto- 
nian. Since the initial systems are in equilibrium, g'^^'> = (g^^^ — 9^^^) tanh Substituting 



these into Eq. ( 88 ) , we get 

I = i J dxdxidu ^ T^{xxi; a;)[tanh ^^y^* {Gif{xxi;uj) - Gl'l\x,Xi;u)} + GI^\xxi;u)] 

(91) 

where 

T(j(xxi;uj) = J dx'dx'iV*{xx')V{xiXi)pa{xxi;uj) (92) 
and the tip spectral function at energy u is given 

pAxxuu) = t[gli\x[,x')-gi'^\x[x')]. (93) 

The exact Green functions in the superconductor can be calculated by integrating out 
the tip using Dyson's equation through the relation 

G{xxi; u) = g{xxi; u) + g{xx2; uj)T,{x2Xs; u)G{xsXi; u) (94) 

where the tip induced self-energy is given by 

\ n ^ ,r(xxi;u;') 
L,{xxi;uj) = I du — (95) 

LO — Ul 



oo 



{xxi;u) = tanh —^^r{xxi;u}). (96) 

The self-energy discussed above is for the normal Green function, which is relevant for 
the tip. However to describe the superconductor, we need to consider the Nambu spinor 
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Green function at complex frequency u. To be consistent with particle-hole symmetry, the 
self-energy in Nambu spinor notation is given by 

( \ 

^Nambu{(^) = \ \ (97) 

\^ -S*(-tJ*) j 

The retarded and advanced self-energy are given at frequencies uj ± i5 respectively. The 
Dyson equations for the full Green function G = {Hsdc ~ ^)^^ in terms of the self-energy 
can be decomposed component-wise^^ as 

G{uj) = {l-g{oo)J:{uj))-'g{uj) (98) 
= (1 + G(^)S(«))(7W(1 + S(^)G(^)) + (99) 

Since the starting system is in equilibrium, g^^^ = tanh.^^[g^^^ — g^^^] and writing (1 + 
SG) = GHq where Hq is the unperturbed Hamiltonian, the latter equation can be reduced 
to 

GiK) = G'(«)S(^)G(^). (100) 
Finally using the equilibrium constraint on the tip, one obtains the expression for the cur- 

m 



ren 



/ = I 



j dxdx.du J2 r<x(a:xi; co) tanh '^^^{Gl^^^xx,; co) - Git\x, x,- uj)}. (101) 



To proceed further we must make assumptions about the Green function for the STM tip. 
For a simple STM tip tunneling to the end of the nanowire, the tunneling can be assumed 
to occur between one point on the tip and the end of the wire. Furthermore, we can assume 
that the local density of states of the tip is a Lorentzian function {{u — /W'^ + 1)^^ where 
W is the band-width about the tip. With these assumptions F is given by 

r 

where /it is the chemical potential of the tip. The corresponding self-energy is given by 

TW 

T.{xxi]u) = (103) 

uj — isgn lm[(jj)W 

where sgn/m(a;) is the sign of the imaginary part of uj. The expression for the current with 
the above spectral density is given by 

I = ^T j rf^X^tanh ^ ^ ^ {^^^^00;^) - C^tHOO;^)}- (104) 
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Using the identity tanh(^) = 4ci;T^^^q Lo'^+(2n+iY-K^T'^ ^"^^ ^^^^ ^^^^ ^^^^ 
analytic in the upper and lower half complex frequency planes respectively, the u integral 
can be replaced by a discrete sum over imaginary Matsubara frequencies as 

(T,n>0 ^ ' 

- TTW J2 ^''''^ ^2f^2^i MG..m ^^t + ^W)}. (105) 
The second term in the expression for / regulates the singularity in the first term for W ~ 



(2n + l)7rT. Using the Dyson equations (Eq. 98), the expression for the Green function G 



can be written in terms of the unperturbed nanowire Green function (?(00; uj) at complex 
frequency uj as 

G,^(00; w) = ^(00; /it+2(2n+l)7rT) [1 - S(00; /i* + i{2n + l)7rT))^(00; /it + i{2n + l)7rT)]"^ |, 

(106) 

The unperturbed nanowire Green g{x, x' = 0; u) at a;' = satisfies the BdG equation 

[{~vdl - - Vg{x) - t^ayd,^-}T^ + V^a^ + Ar^. - uj]g{x, x' = a;uj) = 5{x - a) (107) 

for X > —a with the boundary condition g{x = —a, x' = 0; w) = 0. Here the end of the wire 
has been taken to be at x = —a. Away from the boundary x = —a and the contact point of 
the tip, X = x' = 0, the Green function can be expanded 

g{xx' = 0;uj) = Y^ ^„C^_„e^"^ for x < (108) 

n 

g{xx' = 0; w) = J2 ^nCi,+e""^ for x > (109) 

-Re(2„)<0 

where 

[{-Vzl + t^ayZnjr, + V,a, + At^ - = (110) 
and C are a set of vectors that are determined from boundary conditions. The quadratic 



era- 



eigenvalue problem in Eq. (110) can be reduced to a linear generalized eigenvalue problem 
by defining $ = 2"$ as 

z<^> = -^cry<i> - -[{-fi + V,a,T,} + lATy - Er,]^. (Ill) 

a a 
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The coefficient vectors C are determined numerically by solving the boundary conditions 

^(x = -a,x' = 0;w) = (112) 
- T-J(9^.5((x, x' = 0;a;)|^.=o+ - dx9{x,x' = 0; w)|^=o-] = 1 (113) 
g{x = 0_, x' = 0;u) = g{x = 0+, x' = 0; u). (114) 

The current I{iit = V) is calculated by using the Green function g{x = 0,x' = 0; u) in 



Eqs. (105) and (106). The conductance 4^7 obtained by numerical differentiation of the 



current is shown in Fig. (15). As can be seen from panel (a) of the figure, the Majorana 



fermion mode at the end of the wire (in the non-Abelian phase, i.e. Vz > y + /i^ = 0.5 
meV) gives rise to a peak in the conductance at zero-bias {V = 0). Apart from the zero- 
bias peak, the bulk states in the nanowire also contribute to the STM conductance at 
bias voltages above the bulk energy gap. It is important to emphasize that, unlike p-wave 
superconducting nanowires, the tunneling spectrum of the present system (Fig. ([I5|)) has 
no state (other than the zero-energy Majorana state at the end) below the bulk energy gap 
A. Therefore the effective mini-gap at the ends is order the bulk energy gap (A). As the 
Zeeman splitting is lowered towards the critical value Vz^c = a/A^ + /i^ = 0.5 meV, the bulk 
energy gap closes before re-opening for Vz < ^/'^^^\-^ = 0.5 meV. For Vz in this regime, 
the semiconductor is in a regular s-wave superconducting phase and there is no peak at zero 
bias. Instead there are 2 peaks which are associated with van Hove singularities of the bulk 
quasiparticle spectrum. These peaks become sharper as the Zeeman potential is gradually 
lowered to zero. 

STM spectroscopy can be used to not only verify the existence of the zero energy Majorana 
states localized around the ends but also to study the spatial structure of the zero energy 



wave-functions. Panel (b) of Fig. (15) shows the zero-bias tunneling conductance as a 
function of distance a from one end of the wire. The calculation shows that the conductance 
vanishes at the end of the wire (the Majorana wave-function must vanish at the end of the 
wire to satisfy the physical boundary condition). The zero-bias conductance then rises to 
a peak and oscillates with an envelope that decays away from the end of the wire. The 
decay length becomes longer as the in-plane magnetic field is tuned towards the critical 
value Vz,c = 0.5 meV. No such zero-bias conductance should be observable for Vz < Vz^c- 
Thus, the STM experiment from the semiconducting nanowire provides one of the most 
experimentally feasible probes of Majorana fermions in condensed matter systems. 
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FIG. 15: (Color online) (a) Conductance dl/dV as a function of voltage fit = V- For the plot we 
have taken rj = Ti^ /2m* where m* = 0.04me for InAs, A = 0.5 meV for Nb, /x = 0, a = 0.1 eV-A, 
and T = 100 niK. The different values of Vz in meV are given in the inset. The distance of the 
STM tip from the end of the nanowire has been taken as a = 100 nm. The plots for Vz > 0.5 meV 
(corresponding to Bx ~ 0.5 T for InAs with ginAs ~ 35) are in the topological phase and show a 
zero bias peak while the plots for Vz < 0.5 meV do not. (b) Position dependence of the zero-bias 
conductance in the topological phase. The amplitude of the zero-bias conductance is seen to be 
localized near the end of the wire and corresponds to the localized wave- function of the Majorana 
mode at the end of the wire. 
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The numerical results discussed in the previous paragraphs (Fig. (15), can be under- 
stood analytically in the limit that the STM tunneling strenght is smaller than the thermal 
broadening (F ^ T). In this limit the current can be approximated as 

I (fit) ^^[^(00; /^t + <2ri + l)nT),^] (115) 

(7 n>0 

which can be cast into the more familiar expression for STM current 

/(/i,) ^ r ^ 1 dco'A^^ico') tanh (^ ^'^' (116) 

where the STM spectral function is given by A^a' (w) = Im{gcra' (00; oo)^. The corresponding 
conductivity 

G(/..) = ^ = r E / du'A^A-')sech^ . (117) 

From the previous subsection, it is clear that the topological phase is characterized by a 
single Majorana state localized at the end. This zero energy Majorana mode is expected to 
lead to a zero-bias peak in the tunneling spectrum. To see this we consider the 4x4 Nambu 
Green function near a; = where it is dominated by the zero energy pole as 

giOO;u)^— (118) 

where = {u'[,ui,UpU^). The corresponding spectral function is 

A^{u) ^ u^ul>6iu). (119) 



G{fit) = ^^ = rsech'(^^) (120) 



Using this in Eq. |117[ the contribution of the zero-energy state to the conductance becomes 

^ = fsech^ 
where F = F^^ \ua\'^. 

Therfore STM spectroscopy provides a multi-facted tool to study the properties of the 
topological quantum phase transition in the wire as the magnetic field is tuned from i?^. ~ 
to Bx ~ 0.5 T (which is significantly below the parallel critical field of thin-film NtP^. 
Firstly the STM spectra away from the ends provides information about the induced super- 
conducting gap in the wire, which diminishes as a function of applied magnetic field, goes 
to at the transition and then increases. For magnetic fields above the critical value, a 
zero-bias peak appears in the STM spectrum near the ends of the wire. Finally, STM can 
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also be used to probe the spatial structure of the wave-function of this Majorana mode. 
Thus, the STM experiment from the semiconducting nanowire provides one of the most 
experimentally feasible probes of Majorana fermions in condensed matter systems. 

XII. SUMMARY AND CONCLUSION 

Let us first recapitulate the most important results contained in each section in this paper. 
In Sees. [II, III] we analyze in detail the Hamiltonian and the BdG equations for a spin-orbit 
coupled semiconductor in the proposed heterostructure geometry in the presence of a vortex. 
Here we provide all the mathematical details, left out in Ref. [5T], which are necessary to 
conclusively establish the presence of a non-degenerate Majorana mode at the vortex core. 
In Sec. [IV], we use the same formalism to establish the existence of a Majorana mode 
localized in a "vortex-like" defect in the spin-orbit coupling that can be artificially created 
in a spin-orbit coupled atomic system potentially realizable in an optical lattice. Here we 
find that, contrary to a previous treatment of the same problem, ^ the decay length of 
the zero-energy Majorana mode localized in the spin-orbit vortex is inversely proportional 
to the superconducting gap. Therefore, in the limit of vanishing superconducting gap, the 
Majorana mode will delocalize over the entire system. We then show in Sec. [V] how the 
same formalism can be used to demonstrate the existence of a Majorana fermion mode in 
a vortex on the surface of a 3D topological insulator, even though in this simpler case an 
exact solution of the BdG equations is already available. ^ 

In Sec. [VI], we confirm our approximate analytical calculations, which are indicative 
of the existence of the zero-energy modes in the vortices in the spin-orbit coupled semi- 
conductor, by a full numerical solution of the BdG equations set up on a sphere. ^ Here 
we obtain the full bound-state excitation spectrum for the BdG Hamiltonian relevant for a 
vortex-antivortex pair placed at the poles of a sphere. ^ In the non-Abelian phase of the 
semiconductor film, this calculation produces a pair of lowest energy states whose energy 
eigenvalues approach zero exponentially with the radius of the sphere. This indicates the 
presence of one exact zero-energy state on each vortex in the limit of infinite inter-vortex 
separation. The numerical calculations also show that the excitation gap above the zero 
energy state, the mini-gap, can be made comparable to the induced s-wave gap A in the 
semiconductor film. This surprising result, which was obtained previously for the proximity- 
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induced superconducting state on the surface of a TI,*^ is now extended to the spin-orbit 
coupled semiconductor in this paper. The enhancement of the mini-gap increases the regime 
of temperature in which any non-Abehan quasiparticle is accessible in experiments by many 
folds (from T ~ to T ~ A). In Sec. [VII] we briefly discuss the parameter regime in 
which the non-Abelian topological state is the ground state in the semiconductor film. The 
associated topological quantum phase transition (TQPT), which can be accessed by varying 
any one of the three parameters - Zeeman coupling (K), chemical potential (/i), and the 
proximity-induced s-wave gap (A) - is a transition at which the excitation gap vanishes 
at wavevector k = 0. In Sec. [VIII], we analyze the interplay of the Zeeman coupling and 
the proximity-induced s-wave superconductivity in the presence of spin-orbit coupling. We 
show that, even though the Zeeman coupling can eliminate s-wave superconductivity in the 
absence of spin-orbit coupling, the latter can give rise to a re-entrant non-Abelian super- 
conducting state despite the fact that \Vz\ > |A|. Apart from the zero-energy Majorana 
states in order parameter defects such as vortices, the non-Abelian topological state in the 
semiconductor film is also characterized by gapless Majorana modes at the sample edges. 
In Sec. [IX] , we use the same techniques employed in the earlier sections to demonstrate the 
existence of these edge modes, which turn out to be chiral Majorana modes because of the 
explicit breakdown of the time reversal symmetry. 

In Sec. [X] we study the proximity effects in superconductor - semiconductor - magnetic 
insulator heterostructures starting from a microscopic tight-binding model. The supercon- 
ductor and the magnetic insulator are described at the mean- field level. We determine the 
excitation spectrum of a slab containing a semiconductor thin film sandwiched between an 
s-wave superconductor and a ferromagnetic insulator and identify the dependence of the 
induced gaps on the parameters of the model. In particular, we study the dependence of the 
effective Zeeman splitting and induced superconducting gap on the thickness of the semi- 
conductor film, the applied bias potential and the strength of the coupling at the interfaces. 

Finally, in Sec. [XI] we demonstrate the existence of zero-energy Majorana modes at 
the ends of a one dimensional version of the spin-orbit coupled semiconductor system - a 
semiconducting nanowire. It has been shown that it may be far simpler to experimentally 
realize Majorana fermions in the one dimensional nanowire system because the Zeeman 
splitting can be induced by a modest in-plane magnetic field, obviating the need for a 
proximate magnetic insulatorJ ^° * ^^l We find that the Majorana modes at the two ends of the 
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nanowire can be probed by scanning tunneling microscope experiments. We show by explicit 
calculations that the Majorana modes at the ends of the nanowire give rise to zero-bias 
conductance peaks in the tunneling spectrum at the ends. These peaks disappear on lowering 
the Zeeman coupling so that the system settles into the non-topological superconducting 
state. Furthermore these zero-bias conductance peaks are found to be the only features at 
bias voltages below the induced superconducting gap in the nanowire. We believe that the 
observation of this zero-bias tunneling peak in the semiconductor nanowire is the simplest 
experiment proposed so far to unambiguously detect a Majorana fermion state in a condensed 
matter system. 

We note here that the Majorana fermions, being non-Abelian particles belonging to the 
{31/2)2 conformal field theory (i.e. the so-called "Ising anyon" universality class), cannot 
directly be used for universal fault-tolerant topological quantum computation (TQC).-They 
can serve as topologically protected quantum memory or can be used in quantum compu- 
tation along with supplementary unprotected quantum gates requiring only small amounts 
of error corrections.!^ Since the topological protection for the semiconductor heterostruc- 
tures in our work is very robust, with the energy gap providing the protection being of the 
order of the superconducting gap ( 1-lOK) itself, our proposed system could serve as an 
excellent quantum memory in TQC applications. In a recent development, Bonderson et al. 
have showrP^ that certain dynamic-topology-changing operations, which are feasible in our 
proposed semiconductor heterostructures, dubbed Ising Sandwich Heterostructures (ISH)in 
ref.l^, could allow fully fault-tolerant TQC to be carried out using our proosed systems. 
Thus, in addition to the fundamental significance of the possible non-Abelian quantum order 
and the existence of topological Majorana excitations, these semiconductor-superconductor 
structures may have future application as the basic components of a topological quantum 
computer. 

Our proposed non-Abelian system is possibly one of the simplest to study experimentally 
since it involves neither special materials nor exceptional purity nor the application of a high 
magnetic field. It is encouraging to note that proximity-induced s-wave superconductivity 
has already been realized in a host InAs semiconductor film,'SSIZQ] -^i^ich additionally has a 
sizable spin-orbit coupling. Experimentally, the only new effect that must be introduced 
to the system is a strong enough Zeeman splitting of the spins. We argue how this can be 
achieved also via the proximity effect due to a nearby magnetic insulator. It is important 
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to emphasize that when the spin-orbit couphng is of the Rashba type, we require a Zeeman 
sphtting which is perpendicular to the plane of the film. This is because a Zeeman splitting 
parallel to the film does not produce a gap in the one-electron ban d-structure.'^^'^ Inducing a 
perpendicular splitting by applying a strong perpendicular magnetic field is not convenient, 
because the magnetic field will give rise to unwanted order parameter defects such as vortices. 
It is for this reason that we propose to induce the Zeeman splitting by the exchange proximity 
effect of an adjacent magnetic insulator (we ignore the small coupling of the spins in the 
semiconductor with the actual magnetic field of the magnetic insulator). More recently, it 
has been shown that, when the spin-orbit coupling also has a component which is of the 
Dresselhaus type, the appropriate Zeeman splitting can also be induced by applying an 
in-plane mag netic field.'SSl Note that an in-plane magnetic field does not produce unwanted 
vortex excitations. The one dimensional version of our system - a semiconducting nanowire 
- is also a non-Abelian system in the presence of proximity induced s-wave superconductivity 
and a Zeeman coupling. It is quite exciting that the superconducting proximity effect on 
an InAs nanowire has already been realized in experiments.'^ In this case, the Zeeman 
coupling can be more easily introduced by applying an external magnetic field parallel to 
the length of the wire because such a field does produce a gap in the one-electron band 
structure without producing unwanted excitations in the adjacent superconductor. This 
obviates the need for a nearby magnetic insulator .'^S' In the topological superconducting 
state of the nanowire {i.e., the Zeeman coupling is above a critical value), there are zero- 
energy Majorana states localized around the two ends. Such zero energy states can be 
revealed by zero-bias conductance peaks in STM tunneling experiments at the ends of the 
wire. All other contributions to the conductance occur at energies higher than the bulk gap 
A. There will be no such zero-bias peak when the wire is in the topologically trivial s-wave 
superconducting state {i.e., the Zeeman coupling is below the critical value). Such an STM 
experiment from the semiconducting nanowire will serve as an unambiguous probe for the 
condensed matter manifestation of a Majorana fermion mode. 
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Appendix A: Power series for the Rashba model at r > R 

Even though we found an analytic solution for the null vectors of the matrix for A(r) = 
in the region r < R, we could not find such a solution for A(r) = A > for r > R. In 
a previous section we claimed without explicit proof that the solution to this equation can 
be written in terms of a power-series expansion for p(l/r). Since we are interested in the 
solution at large it! we expect a power-series in to converge. To generate the equation 
for the power-series for p it is convenient to shift to a basis where the 1/r — part of the 
matrix is diagonal. 



Vi-d^r - 4^ + 2-^9, -z') + V,-f, ±A + a{dr + i:-z) 

TA - a{dr - ^ - ^) Vi-d^r + + 2-^9, -z^)-V,-fi 



A + B 
(Al) 



where 



=FA + za 




(A2) 



B = r]{-d^ ~ 1^ + 2-2^^r) + laaydr + aax/2r - ^cr^- (A3) 

If 5" is the matrix of eigenvectors of A then A — SDS~^ where D is a diagonal matrix such 
that D22 — and Du — Trace{A) — rj. Thus the relevant system of differential equations 
becomes 

[SDS + f]{-d^ + - ^<^z + '^zdr) + laaydr + aa:,/2r - ^(Jz] 

4r2 2r2 2r^ \ P^(l/r) 



^-1 ( Pt(l/r) 

1 1 n 

0/ 4r zr zr \ p^_^(l/r) 





(A6) 
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(A8) 



where in the above we note that it is sufficient to provide a series expansion for pi(l/r) since 
S is independent of r. The matrix S may be written explicitly as 

/ Vb^ -c^ + b Vb^-c^-b\ 

where a — rjz^ — /i, b — and c = ±A — az. 

It is easy to see that our equation is satisfied to order by the choice 

PtA^r A ( a^^^(l/r) 

The upper and lower components of the above solutions are different orders in 1/r. Therefore 
we redefine our spinor as 

'pt,2(l/r)\ Pt,i(l/^) 

the above solution at lowest order motivates us to modify our ansatz so that the upper and 
lower component are of the same order in 1/r as below 

m(''''^''A=0 (AlO) 

where the matrix differential operator M is given by 

^ = K ^ I [('' " I + rj{-dl + ^-^S-^a,S + 2zdr) + aS-\iay)Sdr + ^S-^a,S] 



(A9) 




(All) 



1 
r 

The terms in the matrix part of the above equation can be separated into 2 categories. 
Those that preserve the order of a term 1/r" and those that increase the order of a term to 
l/r"+^. The terms that preserve the order are contained within the matrix below 

+ r]{2zdr) + aS-^{iay)Sdr + — 5"V^5] 



J\p. — q2 _ Jj 

+ {2r]z + aS~'^{iay)S) | " " 1 + a- 
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The matrix written explicitly here preserves the order of 1/r" while the rest of the terms 
generate terms of order l/r""^^ or higher. We can check this by applying the above matrix 
to a spinor proportional to 1/r". The resulting spinor is 

- c2 - 6 




+ {2riz + aS~^{iay)S) \ + a 






(A15) 



+ {2r^z + aS-\iay)S) I " H + ^ 1]. (A16) 






From this it is clear that the action of the matrix differential operator M on the spinor 
(Pt; Pi)^~^ is the non-singular matrix that acts on the spinor (p^-, p^) in the above equation. 

This fact can be used to create a procedure for generating the power series for p2{l/r) 
iteratively. To see how this is the case consider the stage after the {n — 1)*'* iteration, where 
the power series has been approximated to 

pt'\Vr) = Y,Pf-. (A17) 

3=0 

To calculate ^2"^ we note that 

(n) ^ (n) 

= Mpi'\l/r) = + Mpi''-'\l/r) = + Mp^^-'\l/r) + o(r-("+i)). (A18) 

Solving the equation to o(r~") yields the iterative relation 

= -Q-i hm r"Mp^""^^(l/r). (A19) 

r— >-oo 

As discussed previously, we already have a starting guess for the solution for n = 1 so that 
the residue is of order l/r^. Using the non-singular matrix we can continue to solve the 
for the higher-order coefficient by inverting the matrix over the residue. The above argument 
provides a procedure for how to construct the power-series for p(l/r) referred to in Eq. 22 
of the main body of the text. 
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